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1. INTRODUCTION 

Applied Research, Inc. is pleased to submit this final report on an added 
scope to Oontract NAS8-34?!57 with NASA/MSFC. This contract has centered on 
tasks which support the Beta Plight Test Experiment which uses an aircraft- 
borne laser Doppler velocimeter to measure the aerosol baokscatter 
coefficient. Tasks reported on under this added scope of work are; 

A. Using the Beta algorithms developed in previous studies the 
contractor will examine the 1982 Beta Flight Test Data for determining the 
Atmospheric Backscatter coefficient. 

B. Install a signal processing/enhancement program from earlier studies 
and provided by the COR on a computer mutually selected by the contractor and 
the COR. 

C. Perform a preliminary analysis to determine the signal processing 
methods which appear to be prime candidates for combining, to yield enhanced 
processing capabilities for low signal-to-noise ratios. 

D. Demonstrate by simulation the combinations obtained in Task C using 
the program from Task B, 
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2. BETA PLIGHT TEST DATA ANALYSIS OP POOR QUAUlV 

Applied Research, Inc. has calculated atmospheric backscatter (3) from 
data taken by NASA in the summer of 1982. See the Applied Research report to 
NASA entitled *'Beta Expeiiment Flight Report, Contract No. NAS8-54337, August 
1982" for a description of the various flights. This data was taken by a 
Laser Doppler Velociraeter (LDV) operating in a single particle 'detection mode. 
The single particle detection mode utilizes a focus position where the 
majority of backscatter detections is singular. A focus of 10 meters was 
determined to be sufficient to encounter single particles the majority of the 
time. 

Applied Research's plan to generate atmospheric backscatter values for 
the various flights was to produce a computer code incorporating a previously 
developed algorithm which required the single particle count data (histograms) 
as input. The output of this computer code is the atraospheric backscatter 
coefficient calculated from single particle data and is called "single 
particle 8." NASA also requested that Applied Research process the volume 
data channel. Therefore, code was generated to also output a volume channel 6 
value. A description of the final computer code will be discussed in the 
following sections and Appendix A contains a complete listing. In order to 
set up the computer code for automatic processing, several computational 
quantities had to be defined. These quantities were: 

1 ) Establish correlation between signal value and channel number 
on the signal processor. 

2) Establish the value that corresponds to "zero Area", 

i.e., determine absolute scale for -§/S. vs. Area. 
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3) Lgff at 10 meters. 

4 ) Waist S/N return at 10 meters. 

Applied Research had previously performed the LDV system calibration 
under Contract No. NAS8-34337* However, operation of the LDV at a focus of 10 
meters required that additional data be taken to completely define the LDV 
system. The additional data was taken by W. Jones of NASA and consisted of 
the data required to calculate at 10 meters and the waist S/N at a 10 

meter focus. Applied Research's "Final Report - Beta Experiment" has a 
description of Lgf^ and how it is calculated for an LDV system. 

Signal Processor Channel Calibration data is required to establish the 
correlation between signal value and signal bin (channel) on the signal 
processor. This relationship was established by W. Jones of NASA and is shown 
in Figure 1. This curve was used for all single particle 8 predictions in 
this report. 

In Applied Research Final Report - Beta Experiment, a plot showing the 

S/N 

relative relationship between Signal/Noise/Cross-section ( ^ ) and Area 

(Figure 2) was presented. Since plot has an arbitrary scale, it was necessery 
to establish an absolute scale for Signal/Noise/Cross-section vs. Area. This 
is the "single particle calibration" which was done using a sandpaper disk at 
a focus of 10 meters. The following section contains the theoretical argument 
for the single particle calibration. 

• Single Particle Calibration 

In terms of the backscattering elements of a rotating disk, the 
signal is S * E where g^ is the gain in the region of the i'*’^ element. 

If elements of same value g^ are grouped together into regions denoted by area 
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elemanta "a"i then S ^ 

differential, then S •» E jiS g da. This can be rewritten aa S ■ p/g(a)de 

e da 

where P “ da bidirectional reflectance of the aurface and ia 

% 

asaumed independent of a. Since g « S/c, Figure 2 can be used to obtain the 
calibration constant C where g «• C*s/<7 • Cg*. C ia some constant which will 
give the proper value to the eq^uation " /Cg'da, where S* is the 

sigrAal return from a sandpaper disk at 10 meters, and g’ is the relativel;^ 
scaled S/N/c' data. Figure 3 shows the calibration data which resulted from ‘ 
this exercise. This data allows a signal/noise/aigroa to be associated with a 
payticular area at a 10 meter focus. 

s 

By using the equation ^|o “ 8(o) to tie.*down the Area curve (Figure 
3) at an Area of 0, a cross-section can be produced for each signal bin on the 
signal processor. This association of ei'oss-seotion with signal bin was 
calculated for each single particle B prediction. 

Algorithm Description 

Vflien all of the required information was completed, Applied Research 
proceeded to set up a computer code which \^ould automate the data predictions. 
The computer prediction code has 3 ma^or sections which are repeated over and 
over until ‘‘.he data ia exhausted. The first section is a subroatine (GETDAT) 
to read the data off disk and translate the proper values to work arrays. The 
second sactioii does the noise calculations for both the volume and single 
particle ehannela of the data. Both calculations are done by using the data 
to "look up''' the noise values in a table, A linear interpolation is done in 
the table to obtain the final values. The third section is the subroutine 
(DATANAD) which is the single particle B inversion routine. 




WlGINfli. page w 
OF POOR quauiy 

18(a)* ijQy j*Q « do ia 8 










The BETA Data is packed into 768 byte records of the SignaX V computer. 
Each record of data must be decoded into the 584 words which are identified in 
Table 1 . The alternating noise samples were used to calculate a "noise 
factor" for each set of data processed by the aigorithm. This was done by 
monitoring the integrated volume P word in the noise set of data and then 
calculating what noise (in dbm) corresponded to this value by using a 
"calibration table look up" with linear interpolation. The noise data blocks 
were further utilized by subtracting off the noise hits in each bin for the 
single particle BETA Data. Care was taken to compensate for bandwidth 
differences between single particle noise and volume channel noise data. 

The Beta Inversion Algorithm requires approximately 10,000 particle hitrt 
to be accurate; therefore, the data had to be combined in most oases in order 
to obtain the 10,000 particle count minimum. When this was done, the time 
that corresponds to the data was also averaged over the data-taking interval. 
Therefore, the output of the algorithm was single particle 3 versus time. 
However, there is not a one-to-one correspondence between input data records 
and output data points due to the combining which occurs when the number of 
single particle hits is low. 

A volume backscatter coefficient was also calculated for the same data as 
for the single particle 3. This data requires integration since the LDV was 
focused at 10 meters where it is likely that most of the time single particles 
occupied the sensitive focal volume of the LDV. A volume 3 is hormally 
defined as the backscatter return from a collection of particles measured at 
an instant. The volume Beta was calculated by using the following equation: 


3 = 


V /V 
s n 


d eff 

where Vg is the integrated volume single, is the integrated volume noise. 
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Table 1. DATA RECORD 


0- TIMH; 

12 

TIMM: 

0 

TIMS: 

0 

TIME : 

10“® 

jTIMHl: 

12 

‘TIMMl: 

0 

iTIMSl: 

0 

<TIMF1: 

0 

8- DATEM: 

6 

DATED: 

1 

pATEY: 

32 

11- 'STEPRl! 

0 

STEPR2 ! 

0 

STEPR3: 

0 

14- VOLBH: 

0 

VOLBL: 

0 

4VOI.BN : 

0 

19- (farms : 

171356 

CLOCKP 

11 . 

FILTER 

2 

20- iTIMEON 

3 

•VCOFRQ 

140. 

iDOFSET 

0 11. 

23- ‘DISCRM 

0 

iIFGAIN 

45. 

25- (fWDID: 

0 

NDSETS 

1 

INTSEC 

1000. 

ADDASV 

0 

ADDAS I 

100. 

30- SCAN: 

0 

TAS ADD! 

: 300 

TASKB: 

300 

THETA: 

500. 

34- VCOOFSi 

! 0 


CURRENT TIME 


TIME AT START OF RECORD 


CURRENT DATE 


STEPPER POSITIONS (NOT PRESENTLY USED) 


VOLUME BETA INTEGRATORS 


CONTROL REGISTER PARAMETERS 
DIGITAL FILTER CLOCK PERIOD 
IF FILTER WIDTH 
VIDEO TIME CONSTANT 
VCO FREQUENCY 
A/D FILTER BIAS 
DISCRIMINATOR THRESHOLD 
IF GAIN 

FORWARD/AFT IDENTIFIER 

# SCANS PER WRITE 

^S PER INTEGRATION PERIOD 
UP IF ADDAS LINK ON 

# MS PER ADDAS OUTPUT 
UP IF SCANNER ENABLED 
TAS FROM ADDAS 

TAS FROM KEYBOARD 

ANGLE THETA IN .01-DEG UNITS 

VCO OFFSET IN 53-KHZ UNITS 


88 IDENT: BLHB 80. 


IDENTIFICATION STRING 


128 

383 DATA 

: ' 9 
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Sj*ie the signal to noise return of a sandpaper disk at 10 meters, P is the 
bidirectional reflectance of sandpaper and is the effective length of the 
LDV while focused at 10 meters. was calculated by using data generated 

by W. Jones at 10 meters which is shown in Figure 4« The corresponding 
is plotted in Figure 5» which shows how this data point fits in with otifitr 
data points calculated by Applied Research. 

The integrated signal, Vg* in the volume channel was calculated by 
subtracting off the average noise signal (v^) for that time interval. The 
signal return from a sandpaper disk at 10 meters at the focal volume waist and 
the bidirectional reflectance of the sandpaper disk were supplied by W. Jones 
of NASA. 

Prediction Results 

A total of 25 flights were processed by Applied Research. Some flights 
had problems with the data which were caused by data transfer, but the large 
majority did not. The prediction algorithm had only two cases where it did 
not attempt to process the data; (l) If the airplane airspeed is less than 
250 knots (such as at takeoff or landing) and ( 2 ) if the last signal bin 
contained more than 5,000 hits. Case No. 2 occurs whenever dense clouds are 
being flown through and since the algorithm predictions are inadequate here, 
the data is ignored. Figures 6 through 27-A show plots of single particle 3 
and volume 3 as a function of time. Figures 28 through 49 show the average 
aerosol backscatter cross-section as a function of time. Figures 50 through 
71 show the atmospheric aerosol number density as a function of time. The 
time on each plot is in minutes, calculated by multiplying the military hour 
by 60 and adding the minutes. Note that Flight 12 is missing completely 
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because none of the data was processable by the singl® particle prediction 
algorithm. Other gaps in the data occur because of the two cases listed 
above. Not© that there are some single particle 6 predictions on some flights 
which get very large. These instances correspond to cloud cases. In general 
those large values of B cannot be considered as accurate since the algorithm 
was not set up to handle high density large particles. The corresponding 
average cross-sections and densities would also be suspect* Figure 27-A shows 
the result of processing with only 1000 single particle hits. These results 
are not significantly different from those for 10,000 hits and indicate that 
shorter flight intervals per output beta point may be possible. Further 
analysis of this question is required. 

The single particle and volume backscatter coefficients differ by a 
constant offset, the origin of which is not presently understood. This error 
may be found in the calibration techniques or processor assumptions. An 
evaluation of the uncertainties in the backscatter coefficient determinations 
would illuminate this problem. 
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5. LDV SIGNAL PROCESSING ANALYSIS AND SIMULATION 


Qti POOR QUALITY 


Applied Research has utilized and modified a signal processing code 
written by the Wave Propagation Laboratory of NOAA. This code was written for 

a Data General machine and utilized on the Sigma V computer at MSEC. 

Conversion and implementation of this code constituted a significant task. A 

listing of this code as it was utilized in this study is given in Appendix B. 

A flow chart of the complete code is given in Figure 72. 

Basically the code is capable of generating a sinusoidal signal with 
noise or a narrow band random process signal with noise, and processing this 
signal according to selected algorithms to give frequency estimates. It is 
the narrow band signal which is appropriate for LDV studies. The code was 
used in these studies to generate the standard deviation of the mean frequency 
estimate as a function of the input signal- to«-noise for a narrow band random 
process, over an ensemble of 1000 different cases for each point. This 
quantity leads directly to the expected velocity error for a given case. 

Use of the Code 

The signal processing program used by Applied Research, Inc. is listed on 
the NASA computer system as ARIDAN and XARIDAN. The ARIDAN is the fortran 
listing of the program, while XARIDAN is the executable source. 

There are two ways of running the program. One is to execute a command 
file labeled DANCMD. This will set the logical unit number ten for terminal 
output. The second way is for the user to set the logical unit ten for a 
terminal output, then start XARIDAN. 

Once the program has been started, data input will be requested. The 
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program will need to know if the user wants sinusoidal or narrow band random 
process synthetic data. The program will then list the available estimators 
and give the user an option to use or not use them. Then depending on which 
estimators are being used, the program will ask for inputs needed to run those 
estimators. 

Inputs that are common to every estimator are: (l) the frequency of 

synthetic data (FREQ); (2) for narrow band random processes, the standard 
deviation of the signal (SIGMA); (5) the signal to noise ratio of data (SNR); 
(4) the number of input data points (NTOT); (5) the length of the transform 
for spectrum (NPTS); ( 6 ) Hanning window for fourier transform (IWNDW); (7) the 
number of estimates (NEST); (S) type of spectrum plot (LOGSCL); and (9) does 
user want a hard copy of input parameters. 

The number of input data points and the length of the transform determine 
the size of the signal generated. If NTOV is less than NPTS the program 
genei'ates a signal of size NPTS. If NTOT is greater than NPTS the program 
Will generate a signal with the size the closest power of two larger than 
NTOT. 

When the Fast Fourier Transform estimator is used the system will ask for 
the FFT averaging time constant. This parameter is used to set the number of 
times a block average will be done. 

When the Adaptive Filter - Weights estimator is used the system will need 
to know the order of the filter (LORD) and also the normalized adaptive 
constant (ALPHA). 

When Adaptive Filter - Predicted estimator is used the system will need 
the same parameters as Adaptive Filter - Weights but will also require some 
others. The program will need to know the number of samples to use in the 
prediction filter. It will also want to know if user wants to cascade the 
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filter with the Kulti-Lag Complex Covariance estimator. This option was added 
by Applied Research. 

Finally, when the Multi-Lag Complex Covarjance estimator is used the 
system will need to know the order of the estimator. 

After inputs, the program will clear the screen and make the plot of the 

first time series and ask user if a hard copy is wanted. 

The program will then proceed to the signal processing and a plot of the 
power spectrum for each estimator will be made. After the processing is done 
the terminal bell will ring twice and a question mark will be displayed. This 
gives the user the time to make a hard copy from the terminal and then press 
return to continue. 

The program then clears display and asks the user if a histogram is 
wanted. Choices of histograms are power, frequency, and width. If user does 

not want a histogram, the user may either restart program without exiting, or 

end processing. 


^0I0R 
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In a previous report (Final Report on Contract No. NAS8-34537, June 1982) 
the application of adaptive filtering to LDV processing has been considered. 

It was seen (Figure 6.5) that an adaptive linear predictive filter (ALPF) 
offers several dB improvement over the pulse pair processor within a large 
range of S/N values for the input signal. However, the theoretical limit 
(Cramer-Rao Bound) for improvement is still beyond this, indicating that 
significant further inprovement is possible. The poly-pulse-pair (PPP) 
processor was know to follow the linear predictive filter result but 
theoretical or simulated results for the ALPF preceding the PPP processor had 
not been obtained. Under the present task, results for the PPP processor (on 
the figure called multi-lag L=4 case), and the ALPF have been separately 
simulated. Then the two processors have been cascaded with the ALPF preceding 
the PPP. As a baseline, the result of an FFT on the input has been calculated 
using the multi-lag L=1 algorithm. 

The processing simulated was taken to be appropriate to the MSFC pulsed 
system, with a sampling rate of 15 MHz. A Doppler frequency of 4*5 MHz was 
assumed, although this is incidental to the result. A narrow band random 
process signal was synthetically generated with a standard deviation of .02 
yielding a velocity standard deviation of .48 m/sec for a total spread of 
about 1 m/sec. A data sequence of 64 points was used, corresponding to a 4 
microsecond pulse with a length of 1200 meters. 

The PPP processor was taken to have four lags, and the ALPF was taken to 
have four weights. The convergence factor, a, was taken to be .1 (See R. J. 
Keeler, NOAA Technical Memorandum ERL WPL-49> P> 15)* 
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Figure 7? shows the result for the standard deviation of the velocity 
error in predicting the mean value for the several methods mentioned above as 
a function of the input signal- to-noise of the signal. One thousand spectral 
estimates were run to obtain each point. The ordinate of this figure should 
be multiplied by .795 to obtain the velocity error in meters per second. The 
results of this analysis are summarized as follows: 

1 . There is probably no significant difference in any of the methods at 
the extreme points +_ 10 dB. 

2. In between the extreme points, there is a region (near 5 dB) where 
the PPP method and the cascaded processor (ALPF followed by ppp) give the same 
result; there is a region (near 0 db) where the PPP is slightly better (about 
.5 db) than the cascaded processor; and there is a region (near -5 db) where 
the cascaded processor is somewhat better (about 1 db) than the PPP processor 
The FFT (multi-lag order = 1 ) is the least desirable in all these cases, but 
serves as a standard of comparison. 

The conclusion from the results is that optimum processing techniques may 
depend upon the signal-to-noise level. In such a case, a more detailed 
mapping of the various regions would be in order. This was not possible under 
the present contract because of the long run times required. However, these 
results indicate that the use of an adaptive filter as a part of the 
processing chain will lead to an improvement over the PPP processor of around 
one db for the pulsed LBV system near -5 db. 

A remaining question which must be addressed concerns the convergence 
time of the adaptive filter. There is no theoretical work available at this 
time on this question for narrow band random processes, to the author's 
knowleuge. These simulations have been guided by results for sinusoidal 
signals* It must yet be established that an adaptive filter can adapt using 
the amount of stationary data available in one or several LDV pulses. 
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:7^.coo c 

_l7t,C00 C.(j TC 30 _ . _ 


l7t. 000 c****** ••• ••••••*••••••••••••••••*••*•••*•••••••*••* 

_J77.C00_C 

1 Ob. COO C 

lOs.COO C *^A)1K'J‘^ f %'T- JKY--YULt/»ALMt ALCCHTt'K 

IbO.COU C 

. I-*! .CCO C . _ . . 

1»»;.C)0 C CFT CCRK^LATinN COrKFlClLNTS FUST 

_ie 3 .tOO ,C 

H4.COO C 


ibt.OCO oOO CALL CCCJKM 
Ibt . 000 C 


_ JK7.C00 C 
iHfc.COO C 

_l^s..cco c 
100 . coo 
1 R 1 , 000 c 
IV/. coo c 
l«5.coo C 


104. coo 
1 vt . co6 



THLV'fTcCCGS’TmSF~c 6FFf'lClrNT£ ’ 
c’all~ CL^ V 

Tll’t~KVsf oV~THr"~VULr/ b ALK F P CALCUL A7ICN C AS bf LOKt 
. SAM* AS THfc OIliER -.AX. fcNTbCPV ESTI.-AlCbS 


cr, TO 30 


J97.C00 C . 
l9t.C00 C 
1R9.C00 C 
ZOO. COO C 
/ 01. coo C 

Z02.000 400 CALL 
Zc3,COO C_ 

Z 04. 000 C 


MAXl^'ll^' KM^LPY- -ADAPT IVF FlLTtP ALGCPITHK 


CALPF 

~ SAVt Trie. PbtUICTtD ' output' IN 'XCUT' 


W07.000 
. 000 

205.000 
21 C .000 
2 n .COO 

1 rv'.'oooT 

2 2 J. coo 

21 4.000 

2lt.000 

2 it.o:o 
217. coo 

> 1 t . COu 

^ ly . 000 


ORIGINAL _PAJ5P rs 
POOR QUALITY 


KKU^ Htkt ON/ 'TIS ThF. SA^t FCK ALL KAX. LNTKOPY 
LSTI-ATOkS 


CCNT I SUE 


CALL NfMKKT 


■join THERtST of Tm£ ESTIMATORS FCR PEAK-PICKING/ 
, I'LOTT I NC^E.TC. 


Cf" Tr ‘0 


adaptive KILTEK--NE1GFTS 


: :-4. coo 
/It. coo 
/It. coo 
I 227.000 
I 2?f'.'CC0 
^iv.coo 
r*: oVco) 

? 1 . C02 

coo 

f/^j.CCO 
■ ?}<.t00 
■'2. •it . OuO 

. 2 3r.'coo~ 

I 2? 7. 000 
‘' 2 ? •? . 000 
coo 

r / < C . 000 
2^ 1 .COO 
2-3/.~000 
r 2 4^. COO 
1.244.000 
2^-t . COO.. 

1 »2 ft . 000 

247. COO 
*24 hVOOO' 

1 .2 - 000 
2 ' u . 000 
•/ ‘ 1 . COO 
. .000 
. coo_ 
1./ -s 4 .'coo“ 
/•>. . coo 
f-.st .oco 

1 2^7.000 
/‘jt.COO 


IF ALPF HA_S_ALREADY BEEN CALLED/ THERE'S NC NEED TO DO 

"■ I T * A C A I N .“. . 

"lF'(IE*^(4))~CC*fO~5f6 

CAL ^_ C ^ Lf^F 

SAVI THE Hrf.DICTkD (jUTfUT IN 'XCUT' 


UO Eu'URlLH A.VALYbIS ON ALIGHTS 

uT'S'll"!' tTaTiT lop d 

_ S!G( IDAA)rACTS(IDAA) 

M)AT = L0PD 
_ MTSF = KFTS 

DC 511 ICA=L0HD^1/4*NPTSF 

SlG(inA) = 0_. 

CALL* SET EH 

FLAC=1. 

CO TC *50 


ADAPTIVE F 1LTEK--PRE.D1CTED 

1 F~A'LPF~HA*r'ALKEADY BEEirCALLED/ THERE 'S" NO ~NEED" TO DC 
IT A r A I N . . . . 


ErtAT = Nr > AT 2 
M'TSF=*FP*TS' 

CALL CALPF 

IF ( IS»^(l2)) g'OTO boo 


DO ECL’PILH ANALYSIS CN EWEDICTEL' OUTPUTS 





c 

t.*\ *COu 620 

000 C 


CALL SUFH 


OKIGiNAL PAGE IS 
or POOR QUALITY 


2c'/*000 C 
4.i *000 C *****’ 
^4. .coo c 

i'ob.COO C 

2 tf *.'000 c 

VI ' I 000. c 

e.COO bO ' 
i"fS0C'L_c 

: jj.coo c 

I'TljCCJ c 

ip.COO c 

-i ?.oco. c 

2^* *1.000 

.•''^.COO 

rr.'cco 

li /,000 c 

1Tb. 000 40 

T<, CoO-C—, 

2-l'.000 C 

: -i._ooo_ c. 

1 Y 

*lsCC0 
V-<.CC0 41 

lyscoo 

^ 1;. OoO 1 3> 

^JsCOg 

y*H.CC'0 c 
iT5.oco.cvrv 
^ -D.COO C 

1 *^ 1 , cuo C 

*■*'^.000 c 

: ^C'^.coo c 

I ' ^fj.coo _.700 . 
^ 1 3.000 

1?7. 000 C 

1 y.b . C 0 0 C 

, J^C0C'_£ 

* 3 u‘0 • 0 0 0 

I ; 01 . coo c 

* 31^.000 C 

t jli.COO c 

' J 34.000 C 

^ -rv. oco 

t ' jI f. coo c 
I jOV.CCO C'*V 
.•fb.OOO C 

jl V. ouo c 

31 C. 000 C 


CtrT‘hr"FTksTinin‘'fKTS OF SPf.ClRUb aVo STCPE THEP 

"C ALV SHFl Kp” ( I ST/ I PL7) 

CHl.Cf' SWITCHES. TO S't_ I E. YOU > AM A PLOT .OF THE__ 

' SP* CTbU*' 

“'lF”(r li-ST) CAlT PLSP£C~ 1ST) 

1 F ( r 1 A C .J Cj_l J CALL PLSPE CllSP — 

FLAGrO. 

CCNTINUE 

hcsTcPt SIGNAL 

'lh' Tf” i’l” I I = i /~i • N p T S E 2 
5 1C( 1 ()1 1 ) = E1GPAS( U'J I ) . 
c n N T I N O £ 

LC 1?- UAl-l/LuPD . 

fcGTS( IUAl)=SAVfc(inAI) 




PO'JllMS FQK COMPLEX CCVAPIANCE 


C A LCULATt COMPLEX COVARIANCE 

LOE 011=1 

CALL CCOkk 


CA"LL”cf)MpLtX COVAOIANCL ESTIMATOR 

*CA fL~CCC) V~ ( 1 ST / 1 RL Z ) 

|»Ll3’T~pbir7F~REG/ A.'ID ”AP HERE 

if' (Ts-^(V ) I’ouYpUTThP'/ ICT/ ' CC'/kSLT( 1 )/ E SL T ( 2 ) / k SLT ( 3 ) 




kOUTINt; FOR MULTI-LAG Cf ^PL.E X _CC V AR I A NCE 




OF POOR QUALITY 


3n.ooo 
31 .'coo 

31 1. 000 
3 147000' 
3 l| ._0^0_0 
i'il .COO* 

317.000 
3ir .000' 
3l| .000 
32C.OOO 
3?,’ .QU O 

|3*:i .003 
3 I j . 0 0 0 
31'4."000' 
000 
v.OOO 
|3 2 7^Cu 0 
{j> .c'oo" 
.CCO 
*3 30.000 
000 
3 3J .000 
- 33^(^0 
3'3j'’.Coj' 
33 .too 
'jjt.OOO 
3 3.7. COO 
000 
CDO_ 
340. COO 
f34|’.005. 

[?4l. COO 
3 4 3, 000 
341*. COO 
3 41 .0)0 
j 4t , COO 
3 4-''.0C0 
|3 4» .000 
3 4^»0OO. 
35J1.C00 
35| .C00_ 
.000 
k 3 ^ 3 . 0 jO_ 
[j^’.OOO 
3‘1..C00 
. 000 
3*;^. 00O_ 
3*^ . COU 
359.000^ 
3tf . 000 
3^1 .000 
jt*;.ooo' 


c 

c' 


C 4LCUL COV A_P AKCt FCh LOkD L ACS 


bOO 


NDAT = M)A12 

'call ccukk 


c 

c 


r.iULhkn SYHMFTRIC CnVAKIANCE FCl 


Lf^DHsLUkOtl 

oOt 1 = 1,LCRDP_ 

S1C(I)=H(I) 

bOb CONTINUE 


M s LlJRL»2 

^'PTS -_LCRD 

CC cflO UM/V2 
blO_ 5 IC( I)=CMHL>rO.O/O.OJ_ 
DC ►'20 Ul/LCRD 
JsM TS^l- 1 


b20 


C 

JC 


f.lC(J)rCnsjC(M 1 ♦! ) ) 

__CU*.T1NUF_ 

*Ub FH ON COV A‘kI an'cf function' 


M>AT = *.nr 
Nkl SFsNPTS 


call tMFi-T 

_GC T(L ^0 



c 

c 

c 

^0 


C'*NT I f.MF. 


A ’iJLtI’ f - l Q Q p 2 


c 

c 

c 

c 

c 

c 


tSD-CF-LOOP 


10 


10 

.t>Q 

c 

c_ 

c 

c_ 

c 

c_ 

c 

c 


is»vl ) = . FALSI . 
F 1 kST = .F ALSE. 


CONTINUE 

CALLALLL 


lilS TCCKAM PLOTTING SEQUENCE 


. 1 - P.L* _1-N_P_A>' I' A T A_ 1 D S F P L J 

2-DC STATISTICS ON THE'm (STATS' 

J-CFT HISTOGRAM FROM THF^ (FHIST> 

4-PLCT ')ATA ANU 'histogram (HPLCT) 


C 

C 

C 

c 


_C.'itC_F_WiyCJiXAPJ.Ab_LF _T0. P0_H1S1.CGJ< AM§ JiVlP 


CALL HPCrFY 
CALL HLLL 


^ftl€4NAL- PA<irt: tS- 

CALL thtJt QP pQQp quality 

CALL HOME yu^>uin 

CALL AS'^MDE 

JIUTPUT. 'i*Q H 5TCCRAMS fi: P : 

CUTr'UT'O-STCP* 

OUTPUT '-l-hESTAKT* 

OUTPUT ' l-fCteEP* 

OUTPUT ' 2-EPEQUtPCV' 

OUTPUT ' 3-bICTH' 

_JNPUT. IMISI 

IS»(t)s.TKUt. 

CALL tPASL 

1K( IhlST.Eu.OGuTtl V9V) 

U ((IH1ST.EC.0).ASU.(1SPC.FQ,1)) GO TO 62 
IF ( IPIST.CI. J) GO TO fc2 
_1F_ tlhlST.LT .Jl)_ go, T0_1 


_ CYCLE jHi?oyCH ALL HJSTCCFJ^«S 


i XC.70 lST = lyt.._. 


C IP UlUN'T JS' THIS L£T1*-AT1CN 

c Gouc TO yo A HISTCGPAm on JT, 

c 

IP(. NOT. ( l5fc(lST7)3 CQ.TO 70.._ 


y • W V W 

2 ^ 3.000 

coo 

396 . coo 


l-.Ki/n J)AIA 

.CALL DSKkU. (ISTy IHJST) 

_ _.3rC.El HiSTOGHAM fpCP THEM 

(.NtTtJ T HE F FT >^AY HAVF A D I Ff £ PF h.T. JiU^rLP 

OF POl.’.TS TO UO A hISTOGRAP OVtP) 


NPl NSsNPTS 

IF ( 1ST.E^.1) NblNS=IJPI3F 

IF (IhlfT.Nt.2) KBINSslOO 
C A LL E H 1 S T (HPA X^rjBl N$y IHI ST) 


^-DC STATISTICS ON THEM 


CALL STATS ( XPEf.Ny STDKV/NRINS) . 


P LOT DATA HIST U.G c A.P . 

r r 7B i = i/Nin . 

blG( 1 )sC''PLX(»GPK(l )/0.0) 

CUNT I SUE 


111 


4 1 b . OC 0 

-;ji.ooo c 

4 1 7.000 _70_ 
41t.OOO 

4r;.ooo 
41,0. ouo 

4 J 1 .C 00 c 

4j^i.000 C 

4 I 3 .COO c 

4.4.000‘"c 

4' SCOO_C 

■4lt,00 0 C 
4*7. CO^. 

4 1' t . C 0 U C 
4 T 9 .COO c 
44 c. coo c*** 

^3l.COO_C 

■47/.*CCO c 

443.000 c 

4‘'4.00D' C 
4>fb. 000 C 
4jt.COO C**’ 
< 3 1,000 C 

'OcVcCO 
-Jy.OOO 
•i'ic.CCO 
<il.COO 

4 T ..000 

•i j_.COO 

■ 4<;-j.ooo 

4jt.C00 

■iJt.CCO 

4n.ooo 

4 f r. coo 
^h.coo_ 

■^'-O.'COO 
4i 1.00 0 ^ 

■ 4i i. coo 

•?*‘3.C00 

■ .coo 

_4r‘_. 00 0 

4i.( ."coo 
4‘'V,COO 

-f-.ooo' 

4 '’•coo _ 

-tO.OOO 

1 » C 0 0 

"4 2.'000 
^cj.ooo r 

" -:;4.'000 C 

4 b.ooo c 

4ot .000 


CHL HFLOT A^/STfiiV/KHlSS/ IHIST^IST) 

CCMP.UE 

CC TC tJ 

i KU‘ 


■ ORIGINAL PAGE 5S 
_.QE_PQOR.QUAUTy 


SOM- J'Jl I*.F INPUT 


"KOUl'lNL TO INPUT’ f POP XEYbOARD I VI fi YTh 1 NC"lT 

UCtSN:.T__K_Nj3.W 


C C ^ C/'S P F. C" / 'J S> C P L ( 7 ) / I S ! ■ C 

CC^^CS mST/* AV,!FST^I«sL'-,LOFUy tLPHA^WCT£{l2P)/N0 

CCV'CS /SP/ ISI»(0:i5) 

CC’“'CN/S1GKL/ S IG< I 024 )/ HSLT( I 024 ), X(JF K ( 1C24 )/ _ , 

1 NPTSP/NPTS, FSTkT/NTOT/ NUAT,PS IC/SS ftf/SNP/ 

21 5YNL AT,GSI_OCNC’lSE/FP_P JySIGVA _ 

“CLPPCN/PLS/ XMJPmTiO)/ ICT/LOGSCL yFIFST/ 1P>FL/LNAPE(4),LPL 

_CrPPCS/LhLS/l AHFL(20yb).._. 

CC^PLtX SIG/kCTS 

LCGiCtL rIP£T,I3#» . . . . 

DIVLNSIUN *'1FILF(4)/IUEAU{256)/lUT(3)/nP(3) 

[1 A T A_ ( ( L A P i L ( I_^J =JL/ 1 0 )y Jf 1 , 8 ) / J_F 0 U ' / J^_R I.E.F. h 

T'" fpA'/'NSFC’'/ 'P m '/4*' '/ 

1 'MAXI ', 'MUM T R 'y.' OP Y-', 'j;BUR'y '.G AL'y '_GCRI '/ 'THM_'^ 

1 ^ ^ 

I'MAXl'y'MUM 'y 'FSTR'y 'UPY-'y '_-YUL'y 'F/XA'y 'LKER'y 

I' ALC'/ 'ORIT 

IJ^AU '^' •^UM 'y 'E sTR',r'OP Y-'/ '-ADA 'y 'PT IV '/ 'E-F I 'y_'LTERJlt. 
' 1 ' ALC'y'RTHV'y 

1 ' AuAP'y 'TIVE 'y ' FlL'y'TER-'y '--HI', 'GHT5',4*' 

I'AUAF'/'TIVE FlL'/'TflR-'/'-PRE'/ 'DICT'/'EU '/3»' 

1 'Cf ‘'P'y 'LEX ', 'Cf!V A ', 'HI AS 'y 'CE '/?*' . '/ 

1 'MULT', ' l-LA’'y 'G CJ'y 'MPLE'y 'X CO ', ' V AR I ', ' A NCE ', 3 * ' 

_F U'J I V A n HJ AJ^( 1 ) y K ( n ) 


GtT SrVf. information' ON THE INPIT'DATA 


CALL AN'^Ont 


467^000 

4^ b.COO C 
4 V,000 C _ 
470.000 C 
'l.OoO C 


2. COO 20 

473,000 

4^'4,000 


_lSHCs 0 

_ _Nb»P-NAPHC*' BAND KAVLEICH PPbCtSS 

DC b£*fcAST SUUSOIDAL UW HhRP OUTPUT OF SYNTHETIC 


OUTPUT 'SINLSCIUALd ) 

INPUT ISVNDAT 

OUTPUT 15YNCAT 


NtHP(2) SYNTHLllC DATA;' 


GIVE THE-. A CHOICE OF tSTI-ATCkS 

CALL i cvVhs (C/tScT) 

JJl'TPn 'Thr FCLLOnING f STI ** A T CP S . API J V * I UfcL E SINGLY O' 
Cv^TPUT 'SI^'ULTASEOUSLY. ' 

_0LTPL7 'IVPLT A l_lf_YQJ *ANT THE iNflCAUP £ ST ImaIOF^. fi 
OUTPUT ' ' 

_r.G JOO I = J^t 

NK17K1O/1CC0) (LAbtL(J/I )/Jsl/l0) 


CDO_1 000.-.f C:-‘'iAT(:0A4J 
COO ICO CLSTINUL 


S'J-( PEPPINT NAME OF ESTIMATOR AND ACCEPT 

__ll 


1 OR 


•.’^ 2.00 7 

CCO 102C 
•]'4.C00 lYY 

*(•'^^ 02 : 

4U.OUO 


i'7.C00 2C0 . CCNTINUL 
^J'b.COO C 


DC 2C0 I = l/t . . 

hPin (10/1C20) (LAbEL(J/ 1 )/ Js 1/ 10) 

C FCFPAT (^0A4/Z) 

Y INPUT 1ST 

CUIPVT 

IE (1ST .E(J. 1) IS«(I) = .TRUE. 


£ST L -- ADAPTIVE (fOUMEP TP.ANSFOFP.). PREI 

IF ESTIMATOR 6 WAS SPECIFIED/ HE NEED TO KNCW HCW 

.'-a_ny ,CUTFUT.J)AIA PClNIi. w^.wjlnuC 

Cr.l FROM THE ADAPTIVE FILTER 


IE (.NTT. ISW(6)) on TO 300 

OUTPUT 'DU YCU_wAM TO ENTER HULTlrLAG fPCH ALPf 1-YES/ 0-J 

INPUT MUAN 

j£ .('iPA!i#ku'.i}i 5 iya 2 i=/ipyE/ 

CUTPUT '*NUPPER OF ALPF PREDICTION SAkPLES TO USE' 

.INPUT M»AT 

OUTPUT NDAT 


IK MAX ENTROPY ESTIMATORS WERE SPECIFIED/ WE ASK 
T_Hi, _I. R_ O RDER. .. 


y ?j^O)_C THEIR OR DER. .. 

*1 4.00J C 

ns COO iOO OUTPUT ' 

^K.UOO LCKU = -1 

4 7 . coo IF (.N"'".(ISw( 2 ).OR.lSw(j))) CO TO 210 

4. ►-nf.fi OUTPUT 'urJtR OF MAXIMUM ENTROPY ESTIRATCRS' 


3 


f 


“ORICTf^TirpAGE IS 
OF POOR niiAi ITW 


Y9.000 
51 0.000 
5*. l.COO 
5/2,000 
4* 3jlC00 
H, 4.000 
525^000 
•Tfc.OOO" 
51 7,000 
520.COO 


•ft). 000 : 
51 7,000 
520.COO 

r s.oc 0__ 

C.OCD 
t ) 1. OOOJ 
‘<2.CC0 I 
y i3i000_< 
T 4.000 I 


! " *;fc. 00 0 ’c 

6S,_C0 0_C 
in. r,nn r 


IKPUT I.UKU 

OUTf'UT LOKi) 

CC_T.r ,32 0 

IF ►*> FNTPOl V M)T SPFClFltD^ PIT IDAPTIVI ISy TH TS' 

WF ASK FUfi THcTk o'kDER (IT'S 1HF SAPE THUG) 

\0 I f“nC)T rr I' sVm )To r7i sirs' ) TorT fs‘i*(6 )T )~ccr f 0 ~ 34 0 

OUTPUT *ORUlK OF AD APTIVF FILTER* 

INPUT LONU 

C uxp UX_L oRO 

CO T'c 330 

i r~'f aTt I V L i r~S P »' cT fTl rTTi'tl. ETf“ SC*'F~bYHrR~NfAT 
P AJL* 

20 I K . NOT . 4J . OR. !S» ' (SJ . OP. IS W { 6 ) )) CC T 0 3_4 0 

30 OUtPl'T 'NCRPALIZED ADAPTIVE CCN's'TANI (ALPHA) ' 

INPUT JkLrHA 

OUTPUT ALPHA 

if' p'a *~e'nTr npY ^ hr^aTa ptiYf " r sT~iVArcP?’Vdf 

S Pjjr_l FJ.5 D AND . 1 F_W UJ^T 1_- L A C CPVP L EX C A V A_R UN C E 

I'S SPEC fFlFD/ 'fri£N’cET NUPBER CE LACS (ORDER) 

40 iE' lLOKD .Et ,’ M~, ANu7'~iy4lTr)~Gb~T^ 

Cr TC 34 2 

Yl OUTPUT 'ORDER OF MULTI-LaG CCPPLtX CPVAPIANCE tSTlPATCR' 

LA I_ L < J K 1) 

OUTPUT LO'kD 

i f~sY n t he' f fc 'dTt'a/ kE neeF“ ro x r. c TY h aIT r eTueYc y / e t c' 

4 2 ^lF~{'lSYrarYTTE'c7o^~GO~T'Fli*5b ' ~ 

PU T P U T_' F_kJ^^L^ CY OF SYVTH R T^I C TA T kj 

input' piYu 

cutput_p;k£u _ , . _ . 

AAILJF _ITAS_NL>RP OATAyAE _NFEP.A 

IJ ( ISVMJA T .FO. 2) CO TO 2323 

343 CUfrUT 'SNK CF SYNTHETIC DATA (Db)' 

1 NplJT SNR 

CUTpUT SNK 

CO TO 3oG _ 

2323 CUTPUT 'STANDARD DEVIATION OF SiCNAL' 

I M'UT _^1 G^A 

OUTPUT s'lG'lA 

CO TiJ 343 


PUT IF ACT U AL DA TA/ WP NEFD SCkp CTHEP PAPAk ETER S.. . 


1 )■ 


114 


ORIGlUAi; 49- 
OF POOR QUALITY 


I 

‘71,000 3b0.. 

‘i'»i.ooo c 

•I 3,000 C _ 
'/4.000 C 
.V^*OQ0_J6.0__ 
a t,ooo 

*l#7,000 _ . 

t7fc, OCO 

^S,ooo 
*1 0.000 
‘.'ll, coo _ 

•i '2.000 
1 3. COD 
‘X.OC'D 
‘•^t.COO Jt7 
i ‘t , 000 

•r«7.C00_ 

‘'iL.'COO C 
• »9, CCO C 
t/0.000 C 
‘.V1.CD0.C ... 
5 2.0U0 369 

?^ooo 

‘<i4.000 3>^l 
^ '*• , CCO 
‘.1 'f .CCO 
‘..97, coo C 
‘^"-‘b.OOO C 

\ '9j.coo_«; 

too. CCO J92 
tOl.OOD , . 

'j *2. 000 393 

'•fj, COO 

'04. CCO 
tj'i‘jOOO j9i_ 
Jjt. 000 

t07. COO 

r'lr.CCO 

d;9.coo 


'395 


39 t_ 

C 

C 

C 

c _ _ 
c 
c 
c 


OUTPUT •‘STI^'^TEU >IDE-FiASD SM< OF DATA (tP)*/$KP 

C^^ HFOPHATION ABOUT .FFT PBOCESSINC ... . , 


>S.T>TM 

IF (f.ST'T .LT. 1) GO TO 3t0 

OUTPUT 'TOTAL .KUMBlK OF INPUT DATA P0INTS'_ 

INPUT NTHT 

OUTPUT NTOT ... ... 

ISPCPLTsO 

n C ’ S » ( u . C 1 ^ W c ^ J • Oil • » s N Q J R £^l 4J . C P .ISh 1 5) . C R a ‘ ( D 


-4 


♦.CR.IS»(U)) JSPCPLT=1 
II- (ISICpLT .Fu, 1) CD TO 367 
CO ‘’■f joy 

CUT! iiT 'LTNCIP CF TP4N5Fri-w FCN 5PICTNUP (PCNER OF 2)' 
INPUT NpTS 

OITPUT NpTB 


IF IT'S APPNCPPIATE TC A$F AfcCUT •I^DLl^»I^C 
(IF IT IS/ DO ID 

’ 1 K'i s ». ( 1 ) 7 0 P*.! S*n ( D 7o r”. 1 S »» ( 6 ') ) "c o'* T O ' 3 9 1 

GO TO J'f2 

OUTPUT 'hlNULtelNC FOR FOUf^ILP TRANSFCKNS (1/0)' 

1 ’.PUT U‘OF _ . _ . _ . . 

OUTPUT UNO* 

and/ L~AST~’bUT "not least..". 


IF ( IS.»( 1 )) CL TO J93 

GO TO 3'^‘i _ _ . 

OUTPUT 'FFT AVERAGING TI^E CONTSTANT' 

_I M'UT NAV 

OUTPUT NAV 

_IF. (I. AV.Li:,().Ji,AVj.l 

OUTPUT 'NU^'tLK OF ESTIMATES 

I NplJT N*" ST . 

OUTPUT r.tST 

_Il- (NEST, GT, 1 000). N‘ST = 1000 

IF ( ISI-CPLT .tU. 1 ) GO TO 395 

_.GC TO. 396 

OUTPUT 'SPICTRUP PLOTS; LENF.AP (0) OR LCC(l)' 

_ INPUT LOGSCL _ 

OUTPUT LOGSCL 

lF(.NbT.Ii:*(t)) V0AT = 64 

OUTPUT ' 1 NPUT WGTS(l)' 

IM'UT NGTS(D • 

output'' i N PUT »iGTS(2')' 

INPUT WCTS(?) __ 

OUTPUT 'INPUT .GTSO)' 

INPUT NGT5(3) _ 

(JUTPUT ' 1 NPUT 'WGTS(4)' 
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fL?3. 

l*?b. 

6 ? 6 . 

n. 

HT 

c?S. 

*30. 

31. 

n/. 

r ?3. 


33. C00_( 
3r.*000 


.'•t J. 

• 7 . 
•w3fc. 
^ IS. 


isrUT hCTS(4) 

OUTHUl'UO YfU h*NT HDCOPY l-YtS/O-NQ* 

IS^'UT HDANS 

IF (Hr ASS.FC*.V)CALL"HDCOPY 
CALL fHkSr 
KLTURN ■* 


original PAC5E4S 
OF POOR QUALITY 


SUtrOUTlNL UlT 


THt CRAM :MTUL171NC ROUTINE 




1 

. c 

no 

cr*'vr 

N/SPI C 

/ 

IS 

PCFL( 

7) 

1 44 

.0 

00 

CCMi'C 

S/IMT 

I 

AL/ 

IRP/ 

IF 

tAt 

. c 

00 

Cl *'*'C 

N Nt.ST 

/ 

r.A V 

/ l^ST 

/I 

U 

. c 

00 ' 


.S/Sl'CN 

I 

/ s 

1G(T0 

24 

47 

. c 

uo 

1 NrTSF 

/ ‘•PI 5^ 

K 

STh 

T/NTU 

* / 

t4t 

.0 

00 

2ISYN j 

AT^CSI 

C 

/GN 

OISF/' 

Fh 

» 4<i 

. 0 

00 

CC^'vC 

»./l*wS/ 


XNU 

>-H( 10 

)/ 

SO 

. 0 

CO 


N/LHLS 

/ 

LAN 

FL(20 

/O 

•-Si 

. 0 

CO 

LCGIC 

AL FIR 

c 

T/L 

PL 



7o 

OJ 

CL'*'PL 

£ X SIC 

/ 

Vgi 

S 


1 

. c 

CO 

ni^EN 

SION D 

A 

TU»< 

S(204 

«) 

1. ‘4 

.0 

00* 

Ol'-h N 

SION I 

D 

AT( 

1024) 



. 0 

00 

EUUI V 

ALF NCt 


(DA 

TUMS( 

1) 

! 

. 0 

00 

DL' 10 

0 1 = 1/ 

f 

CO 




. c 

00 100 

b' c ^ ( 

1 ) = P. 






2e )_/N0 
)/ 


,/LNA^'E 


I'.v. 

6C. 
» M . 


SYMhtTlC UATA (NRk'E) DGES'M MED ANYTHING OCM EOH IT 

N0RRTLT7rf~MULTlPLICAT10N~F’ACTCRS* FOR~SNR 

TOTAL INPUT POWEK (VARIANCE) IS UNITY 


r » . 


SNPKslC.** (5NR/10. ) 

I’NPI SErl . /( 1 .♦SNRF ) 

F5ICsl,-PNOlSt 

CNOI S»=Si;RT(l‘N(JISF/2l) 

CSlGsSuRTO'SIG) 


( 71 . 

f-: 

t7^. . 


GFT.N^HLR OF_TKAIMNG POINTS (rC)__FOR ALPF 

Nf = N'Tin-NUAT 
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ORfGfNAL PAGE 13 
OF'POOR OUSUTY 




f'^ttCOO 

j ?( .coo 
h n,coo 

pt-.COO 

_| 

i ‘40. 000 
.f »l|,C00 
I 1^.000 

.1 t?.000 

ts^. 000 
. . '■’h»CG0_ 
. CuC 
;-7.CD0. 




C' 


.C^T PAMV/TKRS iNTin-FFT 


r-FTSFrM'TS 

IF(NPT5.C*.MCT) CO TO 600 

**_s t L CC U)AT(K TQT)«l.)/ >nC(?. \ M. 


SKT HApAMETEhS FUH PLSHF.C 


► IJ^Sl s.TkUE. 
IPX4-L-0 


^ 1 ^ • C w 0 c 
. ■'V.CCO. 

i^o.cco c 


__XK0f.^(l )sPSlC«NTOT_i .P*mSE 


_ X^CP*'(^)s(S^P^ no>’Q.4 1.) •• 2 


^ '».b_c ^ 0 X V C ) iXAC P PI 2, 


COO 
I »3. CCO 




I »3. CCO 
t H.COO 
j i^J.CCO 
I .COO 
t_w7. f C O 


xNni-p(^)sxscFP(2) 


. fl 


t 


J -n.cco 
CO) 
/CIO. cco 
_y-0i.0o0 

J 02 . coo 
_5-034.COO 
lo'^.ooo 
_ I o?>.ooo. 

1 Ot. 000 
707. coo 
f'Dt.COO 

j?5-a.c.:o 


,/Lnpi) ♦ J./SVPF 
► AC*. ALP^ teLlGHT VtCTJrt 

_P)irXX*«K 

’V*.lli>P(t')sP A’ LCrU 


^ > = (XF*LGnL»)**^ • PSIC 

AlPF GAIN*'} • PSIC . ^ _ 

XSOJ-P(t)spX*NUAT 


XS0F*'(7)sl . 

XKOPM( )sXfiCPP(n 

XNJKP(V)sl. 

XNOK*'( J0) = 1, 


n 


710. OuO 
1 . C00_ 
h 2.C00 

ri3.c:o. 
7H.000 
It.C.C.D 
ilf>.COO 
717. CCO 

I lf. 000 

iC.CCn 
72C.C00 
|:'i.coo_i 
■| . 21. coil ' 
7 2;, CCO. 
r:s.c:o 
I .ir-o 
7/t. coo 


i:LTUPK_ 

£ND 


SIBKOIJTINL CtTSlG (SlGPuS) 


SUBPCUllNK IL GET Trit APkKOPHIATE SIGNAL TO BF PROCESSED 


THtRf AKr THPF!: DIFFERENT OUTPUTS OF 5IGNALS 


3VMPETlC/SP.US7IDAL--*ITli RASDCP NOISE 
(lSYSr»AT=U 


; 3 : 
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ORIGINAL PAQI If 

OF POQR-OUAUZy_ 


.*>?7.000 

I '?e. cuo 

J[.29.000 
730^003 
J •31,_000 

1 »y2.cob" 

_733,000 
|13>,'000' 
|'3*).030 
73t. COO 
f7 3?^COq_ 
■| ly'r . COO 
/3Vj.000 

r ‘ 7^0. CC0‘ 
1A \, COO 
742. OOO' 
743^C00 
J /V4 . O'OO" 
-»74tj^OOO 
”74t . COO 
ri47. COO_ 
"iMt ,000 
_7 4 9^C00 

I 75'C. CCO' 
07 
000 
^7‘3.CC0 
1754 .CCO 
^755. 000_ 
7*^ (.■'coo 
J7b7.000 
475t. COO 
759. COO 
T750. COO 
J7^lJ^C0 

7t'2. coo” 
-7t3. 000 
' 75 4.‘C00‘ 
l76‘ ,C00 

r ut . 000 
7t7jJD00 

T^Tt. coo” 

7t9. COO 


SVMJiLTlC/KbKP--NikRCW-PAKD K4MCK ^bOCESS. 
fclTh 

(ISYNDAT«2) 

KSUlsO z> CENI.KATf NEb DATA 


CC4bC>. HOCK i4lTh SICKAL PASAktfEbS 


CryvON/f^Ji/K^s^SU*-*! 

Cr»'KCf./i-LS/XNnP5{10)y ICT/LCCSCL/ IPXFL/l NA^'K 

CUX^'i-S Kk Sl/NAV/lbST/ 14SD#/!,0PD/ ALFHA/kClS( 1 

c t N__/_s ir.M/ s ir.d 024)/ f h t ( i 4 > / h m ( ic^^ 
“l irp T b k'/ NP T S / A b T P r , N T OT / N y a'T', P S'I G/ S‘KP> / S N h, 

_2 L i> V 5 y A T / C S I C y G N CIS Ky F P_K 0^ S IC.»- A 

Cf'^PLIX sfc'/bCfS/CSCAL 

_l».TkC^b haYli _ 

L'lMf NSIUS DA1 Om 5(2O40)/ IDAT( 1024 )yftAM 1) 

KUl V ALk SCJ JJj A TUMSd )^ SIC(] )) 

log"kVl*f fbSl* 

DATA ‘.Sktryc/ _ 

l'F(?.Skt D. bt. OCALL RNUU(PtS, SStLD) 

c’o~T cTT pf pjFh iTt l ~p o u t T h 

ISsISVM^ATTi 

CC T C ( 2 0.0 Oy 2 p pPy 30 00 J. .1 S 


4),LPL 

2H)/b0 

)/ 


J 7 ^lJ^C 0 

7 t” 2 . Cbo” 
- 7 t 3. 000 
' 75 4 .‘C 00 ‘ 
l 76 k .COO 

r ut . 000 
7 b 7 _dD 00 

TFt. cbo” 

7 t 9 . COO 

F 77 0 . 000 “ 
77 1 .000 
" 772 . COO 

{ 773.000 
■ 77 ‘ 4 '.'COO‘ 
* 775,000 
r-77t . coo' 
I 777 . COO 
77 b.'COO 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

2C0C 


SyNThiTIC DATA 

(SINUSOIDAL) 

CFJ. SUE_WAVE FROV DISK __ 

(IMT HAS PLACED COSINE IK PEAL PART/ SINE IN 
IPACIN AHY PART) 

SyslHmC_DA_T_A_(SI».US01DAL) NkE_DS A SINE NAVE IN *SI KE 

(PE'a'Ls'COSINE/ IRACINARYsSINE ) 

~CUTPU7 

_ IF (k >j)4C 1/ 402/ 402 

F‘PEu = kRL'.”»l . ' 

Ue6, ^b3ld50^*k'KEU 

'D0*4 1 o' Uf/NkTSF' ' 

TPLTA sK-FLOAT (1-1) 

“S1G( I ) = C‘'FLX(CflS(THETA )/SIN(THh7A )) 




*). OOD 


).|C.O.O)OUT?UT l/THFT* . 

crt.T isUf 

TMls^.O . . 

F I^STs.TKUe • 

S ! C» * \-j 1 ♦ 

CM* I . 

► s 7 ■ 

ADO VaSDOM* KO!SE TC EACh PAH‘ 


ORIGfNAL PAGE IS 
OF POOR QUALITY 


k 3 C 0 


N’OTisM TSF *;■ 

i»c :^ 3 oo isi/Mm ^/2 

i Au*'r.( I ) = L Aiu''‘‘.( I )*csK»rpAN(LU‘-)*ci.cirr 
L*Tl^^i,(l*l)sLATU»'S(l»l)M;r.IC*LFANCC)UF)*CSOISF 
31 Cf'*K = l» 

CC T*" AjOC 


C 

£ • • • ■ I 

c 

c 

c 

c _ _ 

c 

c 

c 

c 

c 


Ci 7 . 000 .c_ 

8 lb. coo C 

(T v,ooo _ 
A 0.000 c 

c^l.OOO C 

tj ’^»C00 c 
(i 3. coo 


Syi.THrTIC L’ATA 

( ^^A>. hC. Dj!,A'nOwjknCtS§3 


CtNrhATE SICSAL Sl'iCTP'Jy AND SIGNAL FO*<ER 
USING FREQUENCY AT MIDPOINT (.b) 

PJO. . 

Csl. /(SIGmA-S(.KT(2.*P1 )) 

SlG2s;;,*SlCFA*5igKA . .. 

PSlGsO. 

DC jilO I = 1/NPTSF . . „ 

LFPM. = MnAT( I-l )/FLOAUNPTSF)-.5 

>SLT(i; = Q... 

IF (AFS(r)KPFL).LT.{5.*SIC‘'A)) RSLT(l)sC*EXP(-(DFPEfDFPFu)/SlCi) 
P£IC^f-SlC♦^SLT(I) _ 

CAUL'ULATL.NOISE PONLR 


.F.nOI 5 L=PSlGi§NKF._ 


N0K*'.AM7.L TO UNIT TOTAL PONEP 


^ I 

X'- 


HSJ-sPSI>>PNClGF 

■pSIGrPSIC/DSUP 

Fr.. 5 !SFifNiJl*IX.PSU« 

= l.C 


'V. COO 
.^C . COO 


RANLCMIZt TO EXPUNFNTIAL PISTHPUTION' 


bCAsl ./FLCAT(NPTSF) 


32 ;o 


33/1 


. 000 

I f'SV. 000 

. CoO 
‘>7.000 
. 000 

1 S‘<. coo 
^ c.'cco 

tt 1 . OOO 
tol'.OOO 

Jc3. coo 

b t 4 • 0 (' 0 

1 0 ^^. 000 _ 
tC, oou' 

»>7.C00 

J r. b . 0 0 0 
. COO 
70. OUO 


in 3. 000 
'b74. 000 


3230 


3233 


3234 


40C0 


^ 

l»n 322 0 I« 1 ,M'TSF _ OF POOR QUALITY . _ . 

NSUlsl 

C»LL KSl»U(l- >S/Nsn- 0) 

PSL1 ( 1 )s*ALCC(»4AN( 1 ))*SCa*(PS 1C*KSL1 ( l)4fl.01Sn 

G *:'nt ita t i> i: t h“'sp f. c t k u o*^ n'v e b s e ' f f t 

ShlUlNC SPFCTU4 nWH so that CkMEB FkEOUENCy IS 
OVtB SPFCIFIEO FHFuUtSCY (FBEt) aND 'TAILS' OF 
THE SPECTPUM art. CYCLIC 

CCf.TlNUF 

T^E*«.2sF^*^w 

U (FHhw.LT.O.) FIE(.2sF RE uM. _ 

f.Lsm«(M-i(.2*NPTSF)-NPTSF/2*NrTSF*2' ' 

__ 

tC 3 2 30 Isl^M’TSF” 

ILsbf I ( ( IrSL )/M'TSF)*l 

"l \U* I s/* ID ■ * " 

|SDXbslf«JXl-l _ ___ 

■xr.sSiRT(KSL1 (f) )•. I 



CALL H.DU(MN^XSFtD) 

AtGU-s^.'Pl'BASd ) 

D AT tJ S ( r.D X F )Vx N • C'US ("a kCU y ) 

DATU*'r.iisbxi) = xs-sis(AKL:u'') 

AVCrtVCytnSF 

itUTFUl AVG 

GfMFATF CUHPLt X T I SJ bJ_tS_ f PC P_ I \VE KJf FF T 

► sALCG(FLOAl(NPTSF))/ALiJC(2.)^.5 _ 

OStls-1. 

CALL FH (‘'/CN-l) _ _ _ 

SCALsKLUAl (KPTSF ) 

CSC ALsCVHl X^SCALy 0. ) 

S-J^-xYsC. 

_nr j 23S XPANsd2*M'TSF _ _ _ 

blG( SDAN) sCSCAL*SIG(NDAN) 

SUMXlrCAbS(SIC(NDAS'))**2^SUVXl . 

K vrs5j^T(SUyXl/(KPTGK*2)) 

SCALsl/kVS 

C'S C A lVc y P L X ('sc A L , 0 . ) 

SU'^XIrO. 

10 J234 NDANs iV2*‘nPTSF' 

_ _ SI'UNOAX ) = C‘•.CAL•S1''.(^^.A^) _ _ 

Sll-X I iSlbX 1 ♦CAhS(S1 G( N’UAN) )**2 

h' S = SLH (SUMXl/( f.?TSF_*2)J 

NnSr=64* 

PrTUKS _ 

LM) 


original PAGE IS 
OF POOR-^iUAUTY — 


I I 

,> ? t • I 

•oth « 

e*'9. 

J^* 9 J . 

e<»2. 

1 ^ 94 . 
_ ? 9 t . 
T - 9 t * 


7- , 

-I • 


000 
CDD 
COO 
000 
COO C 

cob 
coo 
000 c 
C03. C 
000 c 

(00 

cob 
coo 
ecu c 
(CO c 
coo c 
Q&0-_ 
coo 
coo . 
cco 
coo c 
coo 
coo_ . 
000 
COO.C 
coo c 
000 c 
000 
000 ... 
coo 
cco.„ 
000 


cco 

CO) 

000 ^0 
(00„_ 
coo 

coo 

000 


CC»'^l 

C( ( 

I KHh *?: 
/ISYM 

rc*-* 

Cl -' i ! 

gfj*'i 

U C I( 
1 _ J ’• 

XK ix 

y 

Y'li 

; I Li . 

rc i< 


*5liCFCUl I Kt’ 1 UG'hK rSlCP4K/SCL)'” 

'c C ^ U N / • ’ L S / X N c *-'h u o') /Tc L 0 C S C l 7 n R S T Tn * F L 
CC^^^^rs/^^K.SL/f ]C(1024)^RS1.T( ]Ci><)^>>PRX(lO;4: 
1 KHr'TSK/ NHTS/KblRT, NTOT^ ND AT/ PS IC/ SNRF/SAR/ 

nSYNDAT/CSICyCsgiSt^FRf.ii/SjCHA 

CC»'^rf. /$•/ lS*(0:i*>) 

lUN 1 


. FIP 

0 . 

► LOA 
-9CL 
:GCL 
T. IS 
'LL t 
iC=l 


1S*(0: i*>) 
/►-► S / 5 UMX ! _ 
LT(3)/ITH(J) 

C ^ hf.lS 

*T/ ISR 

1 (NPTS) 


>iCU 

R ASL 

/ i 


RETURN 


. IF K =1 PLOT REAL LATA. IF 1C.* 2 PLOT COPPLfX TATA 

._itLL..T*l».:APitO>9eo^ll0^7(.0) 

IF (IC.i-v..*’) CALL T.«lNLJ(EO/9^.0/fcO/ 410) 

CALL Pi^lRL'J (AMN^F.P.AX/ VP r*/ Y^'AX) _ 

.. . _ ['RAX AXIS AND TIC fIARFS ... 

_i A L L ..Pi ! V L A .. L ' J 1 1 N J _ 

C«LL I'RA^A (XPIN/Y'^AX) 

_CALL PCVLA (C./O.) 

CALL DHA«A OPAX/O.) 

'"call PLVE*A~('xprN7YMrN~) ” ' 

_.C<LL_St.^L"C iJyylYJ 

CALL Dx^AdS (IXtjO/IY) 

uhaVti'cs ~ 


NY s 
.YV .r 
NVPl 

Y INC 
UC 2 
CALL 
CALL 

_ c^ll 

YV = 
YV = 

Y INC 
I'L 3 


10 

_ 0 . . 

= N 

= ( 
0 J = 
_ KPV 
' SF.r. 

Y V 

0. 

= ( 

0 

PUV 


Y ♦ 1 

YPAX C.)/N'Y __ 
1/NYH 

i,A iC.yYY) 

LUC (IX/IY) 
Ai ‘ S_(I XO _^ JYJ _ 

♦ Vise 

0.~-'"i‘MN)~/NY 
r A ( C . / Y Y ) 
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•»H7.0QO 

COO 

VHl .000 

V9X.OOO 
VVt.00'> 
S'/. coo 

‘^‘S^J.OOO 

90^.000 

9<jT. coo 
994.. 000 

997.000 
99T.OOO 
991 . COO 
000*000 
OOf .COO 
OCj .000 
C(Tj._COO 

o'cy'.coo' 
0'-[ .000 
ool. coo 
007, coo 
^ ocT.coo 
: oo , coo 
^ 01 c. coo 
0 1 ] .coo 
ClJ .coo 

013.000 
ojy.coo 
0 jJ .000, 
Olt.COO 
OL7.000 

o:| .coo 
or., coo 
020.000 
0 n J 0 0 0 
”c J. . coo“ 

02J.OOO , 
CY'«COO 
f .o;^ .COO. 
* 026.000 
r iQy»,coo 
>o;j *.000 
I 10?9.000. 
^ lO^P.OOO 
i lo:? .000 
f lO'h, 000 
I 10 33. COO 

f 10 ri.coo* 

IC1..C00 
' 1036.000 
1017. cun 
,lO 


^ - ORl€lMAL" fm'f- IS - 

^ OF POOR QUALITY 

cc'-‘'r s7 s iTs lTTi c < i o 2lT7« s l tTi ?2'<T)7wFp 7 rro'iT)'^ 

1 KPTCK/f.^'TS/KSTKT/NTJT/NDAT/HSIC/SKkF/SNh/ 

21 SYNLAT/CSIC^CNOISt /Fhfig^ UICHA 
_ _ CC“‘'CS NcS1yKAV,IFST/I«*.bk,L0SD/ALPHA^liCTS(l26)/N0 
cc»'^CN/sw/iSi>(0: 15) 

LlJCICAL FlkSTylSW 

CLkPLtX SIC/fcCTS 
C 

C KlGLkE UUT 2 **MsM'TS 

C _ _ _ 

k=aCC(FLUAl(NPTSF))/AL JG(2. )♦.* 

C _ 

'C " “~U0 *1 i»ANT *1900*1 SC ((UKT**!/ ‘yes')" 

C_ 

li«liC=l*NL-* 

Jt- ( ISK('i)) CCTO 1091 . 

GOTO 1092 

_J09l U0C = 0 

10»2 IF (UOO.NF.l) CO TO fod' 

C 

C YcS 

C _ __ 

C FICL'F-t OUT isCWE'^FN' ^ Ch THETA EETkFtN EACH 

C ELtHEj^T 

^ - - - . - 

WFA(;=C2.'3.Hl59265359)/r’LOATa;DATM) . ...... 

C 

C PROCESS ENTIRE ARRAY THROUGH HAN9TNC WINDCk 

C (INCREMENTING ARGUMENT TO CCSINF (THETA) BY 

C *FAC EACH TIME) 

f. 

ThETAsO. __ 

CC 10 I=1/NCAT 

TI-ETA = THcTA4»F AC 

SIC( 1)=S1G(I)'.5*(1.-C0S(THETA)) 

C 1F(S1C( I ) .F(..0.0) UUTPUT S I G( I )^ THET A , *F AC 

10 CONTINUE 

C __ 

C KILL CUT RITH ZEROES IF NDVtVnPTS 

C _ _ 

100 ~ Nr,ATl = NDATM 

_ ^nr no i sndat n npt9F 

110 S1G( I)s(0./C.) 

c 

C ’ ‘ U0"FFT UN SIGNAL 

C 



H. 

J 




( 03V. 000 
040. 000 
0^1.000 
104?. 000 

ro4 J. con 
~io4 4,000’ 
104t>. 000 

r C4fc.C00 
C47.000 
I04h,000 

in?. 000 
.1103. con 

Tl 04 , 000 

t l v, . coo 

iGc. 000 

_.lCOj^COO 

‘"I'lOf-.COO* 

.Ins, 000 

41 10. 000 
illl. CoO 

1 ! ;> . coo 

.»ii.OO0. 

J 1 I 4. coo 
Y 000 
lit. 000 
.ill7. 000 

11 It 4 coo 

_[i IS. 005 

A?o. 000 
ir;i. coo. 
;T 2?. 000 
4 23. coo 
1124. Ouo 

J 2Sj_0Cp_ 
2t . OuO 
.1127. coo. 
;*t’ . Ouo 

.i 2^. coo 

A 30.000 
,U3UC0')._ 
2' 3 2. Ouu 
n33.coo 
34 . coo 
1} 3t),000 
K3fc. 000 
1 1 r/_j^c'>o 

irut-. coo" 
11.3V I 00 0 . 
J 14C. coo 
:j'4j , oc 0 
1 4^. coo 


QRIGINAb PAGE'IS 

OF. POOR QUALITY 


c 

c 

c 


CALL (H^rKF) 

PLACl_ P_n.«.FH SHF.OTPUM INTO, kLSULT 

nr. 2C0 1 s 1 / nF'T'*.k_ 

NSLT( I )sCAB'S(SlC(OyV/rLCAT(NDAT) 


c 

'200 


c 

c 

c 


c 

c 

c~ 

c 

c 


_ IMN1»LT( I ).tt'.0.0_)_0_UTI*UT I / S 1C( I NDAT/NSLT ( 1 ) 
CPNTlNi'F 

_ TUkN 

tM» “ 


SL..)Cf|lj1 1 St C ALl-F _ 

ALHl (AUAr-TlVt LiNbAh PK kM CT I ON* H LTE N ) 

ASPAV * I n_TO _UHiATb CnfFFICUNT VFCTCk C ACCCPDINC TO THE 
L^3 (SUSY CPAUtST) ADAPTIVE ALCUklYhM 
FJLTIK IS FILLtU bb KQkS ADAPTION bEClNS 


C 

C',-**0n/LA/FLACQ . 

CC-*'cr./S ICNL/SI0( 1 024),PSLT( 1 024) / NCPN ( 1024)/ 

I K I TS b / N T £ y A ^ U T^ N TjriT^.N H A T y I* f, I Cy S N KF/ S N 

2i£'vM'AT/<:bif/Csui3F/FPt 1 / si'cp'a* 

CC-^.S Nr ST/MV/ 1^ STy U’lD.%/1 OwD^fcLPHA/»CTS( WP)/NO 

A SIC /*'.TS / XUbTT,tr</t'POK 

X-L' = ALPHA/h LPATCLCPO) 

CC 1 1 - 1 / I.UkU ’ ' * ' 

>.r.TS ( 1 ) = 0._ _ 

LCKUP = LUK’l»tl 


UJ TKAIMNC 

DC 100 KsLDMJpy’ 44d' 

A':o-i.T=(o,y p^; 


10 

c 

c 

c._ 


FlLTfk TU CET FTH OUTPUT 

PC JO LM/LCPP 

M( = K-L 

JL' _y T T = X i I U JUl ♦ W C X S i IJl* JJ.C i K K 5 


CL'NT 1 NUt 
'ad apt" r^MCHTS 


1 


M- s D I C ( r: ) - X L U T T 
t POh = *•> • X'^U 


20 

100 


Lfj 20 L = 1/LCPC 

. ,FA = A- L . . 

N C T £ ( I. ) = «GT 5 ( L ) ♦ L kOK • C fiS J C ( .£ f C ( X K ) ) 

. CCST I M)K _ _ 

CCNT P.Ub 


{ 


A 

- - ! 
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43.000 10: 

44.000 C 

4*5.000 C 

4t. OOO' C 
4V.0G0 

Ve.coo 

4V. COO_C 

5C.C00 C 

51 . 000 c 

52.000 

53.0 00 

54 . COO 

55.000 /lO 
5t . COO* C 
57. 000_C_ 
5b, Coo C 

^5^000 

'60.000 

6l_. CCJ 

62.000 
63^00 

64.000 C 
65.CC0 2^0 

Vt . 000 c 

^7 . 000 . 
be. 000 200 
59. COO C__ 

Yo.'ooo 

71.000 


Korsf;o» j 




T H A m sc O VFk/ SA U OTFUT 

L'C_20£ Ks4 4S/5l2 

"AC'JTTs(b.,0.) 

*H L T E K T 0“C E Vf iToU fp U T 

*C'C" 2 1*0 "L*='i7 L CT i) ~ 

Kt'sff-L 

"X LO'T T = > 'JO T T ♦ h C T S ( L ) • S I C ( X K ) 

CCNTIUJF 

a: APT_»f. ich_is_ 

E C ( K )-X CUTT 

IXOP si X 'XX-U 

rc 220 Lsi^lCKD 

Kl sK-l. 

.CTS(L) = «i(;T5{LWr.kGfi«C0SJC(SIC(xr)) 

OUTPUT lAcTS(L) 

CPNT INUK 

"l " D s f. vV ♦ c a •< 6 ( f k ) 

riC(K-4-'tJJ_=_X0^TT 

CUSTISU5 

11 ( S T (’T ,_N C . KCJ _E^lirM S./ ( NTO T-S C) 
PrTUhN 
LSD 





ilC. OUO 

i Jll.coo. 


f.U(iWUUUNf CCCJHK 


"CliMJLkX CUKKI.ATYuN hOu'nNE 


t» h; . 

LUHC 

'ntoY 


• wwrr 

.-SICNAl . .. . . . . 

Ufvutk 

ou Ti'Ui_ rn«;F F l cifnts 

NUMHLK Qh 'lOTAL DATA FOUTS IK SIC 


.QRlQJWAL’-P^lfHF rs 

^ POOR QUALITY 


_» Jl.v v. 

( 3 J • 00 ' 
■J^V.OOi 
J30.00 
I 331 t 00 i 
I COi 
jjlUOCK 


cf - -cs/r-zis^coMS) 

__C'>*'CN .N£.5T^_**_AV^iLr>T^l-VDwyL0FDyALPHA^>CTr(12U^^fi- 

CLK ‘CN/0 loNL/ S 1C( 1 0*< )/ «yi>LT( 1 024 )/ bijK «^( 1024)/ 

lNFTSF/Nl'TS/^S7/^T/NTOT/Ni)AT/l’i;iC/SNHF/SNh/ 

21 bVSOAT/ CSlt/CSOlSbl/F^L </ S IC-A 

. CC*-KL^x Sir.^KU2V)/C5UH^k»f.T$ . _ 

fcOUlVALtNCt (K( 1 )/l^SLT(l )) 

LOClCAL_Jih 

ILUftlt = Li)KU 

If's’i'c AKWaV fs^COM'Ml KfTFI<nH Y'lPf” US t~THt**F irst^noat pt*s of THF AF 

~ IF (.NUirdSii’dY)))’ CiJTO lO^j ' * 

_Kr;.= i 

NO Js Mj AT 

. .CUTi) 12Vo 


C 

c 

c.__ 

c 


> c 
) c 
» cc 
>. cc__ 
) cc 

I 


) :m 7. coo _ lObO 

J 4 t .000 

.MSCOO 

IJbO. OUO 

tjbl . coo 200 _ 
1 3 S 2 . coo' 100 

fj 53 ,COO 

»j 54.000 


IF Sir ARRAY IS KPT-.CC^^IKC fRC** AL?f .USf 7»'t LASI.KDAI FTS 

_ sncs 1 

NUCs'.bAT 

10 100 L = 0 / 1 LOJ-P _ 

CbUM.sC'>»*LA(C./0.) 

NVL = K*>fl“L 

DO 200 KsMJL/NHL 

_csu’< = cr.ijM»rcNjGf£i CCFJ )• Ji c<k»l) 

R (L» 1 ) = CSUM/MOT 

_R^ TURf; 

tRD — ~ ■ 




403 . 

I 

I 06 . 

40 t. 

I Ofc. 

40 S. 

« 10 . 

\ n. 

lii. 
[h. 
t . 

AU . 

ri7. 

lib. 

41V. 

1 :c. 
21 . 
22 . 

1 23 . 

]"• 

42^. 

<?V. 

] 3 C. 

43 U 

432 . 

1 33 . 
34 . 
3 S 
lit. 

4 Jt . 
i 3 v, 
j 4 C. 
'41. 
441 . 
|4 3t 
1.44. 
44 t . 

f4V. 

|so. 

*'jl» 
4‘.2 . 

[V:-. 

1 . 64 . 


00 0 
000 
000 
coo 
coo, 
000 
coo 
coo 
coo 
cco 
coo 
coo 
coo 
coo 
cco 
coo 
coo 
coo* 
000 
cco 
oco 
coo 
000 
coo 

C"0 

ooo 

coo 

coo 

coo 

cco 

coo 

coo 

coo 

coo 

coo 

coo 

coo 

cou 

coo 

000 

000 

coo 

coo 

coo 

C J 0 
CoO 

coo. 

000 

coo 

cco 

CO'I 

OUO 


c 

c 

c 

C- 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C*' 

c 

r 


__ _ SII-:»KUU 1 I Nl- .(IST^J) 


ORIGINAL PAGE 13 
OF POOR QUALITY 


RUUlIKf TO CALCULATE FIPST CIPCOUB VOM£KT CPEAN 
FHilUENCy) AND HALF PCiiLk klUTH CVER INTIkE SPECTKUk 
AND TH.^N PFCALCULATf PPNFR CVIP C^LV SPECTRUk 
bANCNlOTH ABOUT BEAN FkEC. 


kSLl AkkAY kUH SI'F.CTPUH 

1 ST NUMbtk JK ESTlkATOk 

.) NU«ttk OF estimation 

NPTS NUMBIk OF POINTS IS RSLl 


C _ 

r 

_ CLV-O*./PFS/PFAKS( 3 /l 0 / 4 ,e) 

CCHM JK /SI GVL /SIGH 024 ), kSLT( 1024 ), kOPM 102 4 )/ 

IMTSk ,NkTS, kSTk T,NTJ 3 Ty NU » T, PS I C, S Nk I , SNh 

Cf MM - f. / kl's/ * N C R M ( 1 0 ) , I C T, L 0 C*SCL , F fk ST, I P tVl 
Ci.MfLIX SIC 
C 

C FIND FIFST circular kJHiNT AND NIDTH 

C 


L k =.t , .-5 il Wf L OATi NflSJD 

SUHCsO. 0 

5t‘MSsO.O 

P. I.FR s 0.0 

FK.MisO.0_ 

kim H r 0.0 

^CC.JC0 IsL#M'T§F 

Slsl-SLT(I) 

wspi L».*F LUAT ( I- 1) 

SUMCsSDMC^51*C0S(k) 

SUM£ = 5UMS^Sl'SIN(k) _ _ 

PC^LR = PU»EP ♦ SI 

Fkukl_l F>.4F 

NIDTH = kIDlH ♦ SI*N*N 

J00_ CCNTINUK _ _ _ 

FF«.Mi=H^..Ml/pr*ER/6. 26 31V 
Fh.MjL' Art TAN2 (SUMS, SUMC)/t. 28319 
C LUTPUT J,FPFtPA 

N F Pu =_F P \ u S A_^N PX5f-tOi.^ 

C 

C JIM'I )**2_= J>9, 4784 

C 

k IPTHsw IDTil/F 0.»F P/39.4734 - FPGM1*FPH(.1 

•iDTH=Sw^T(k ILTH) 
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I 

I 


I 4St.C00 
M( t.COO' 
4J l.CtO 
^Sb.OOO 
<fS^COO 
4| O.OOO' 
l.CuO 
) 4i/.000 
^ 3. COD 
1^4, COO 
I . CwO 
4jt‘:COO 
0*^0 
4c»-.c:o 


4 Ts.o:o_^( 
cl C. COD 
147 l.C 0 D_C 

i^rrcoo c 
:<] j.ocD c 

1O4.C0D 

1451.000 

njt.coo* 

' 14/7. COO 
i*-;7c. coo ■ 
f :4|<f. c. 0 

■ j^ic.cco c 
' 24 <il.c :0 c 
^ o4l:'.coo c 
I 34fij.C00 
( l‘4c4.*C00"C~ 
r 3 <|t.CC 0 C' 
’ l4|t.C00 C 
: :-3-f7.CC0 C 
. 1-|*-. 000 C 

1 •; J V . C 0 o_ c 
:4 ^o 7 coo c 
r 1 4 ^ 1 . C C 0 C 

i :<} ^.coo c 
f 2 V^.O. coo C 

24^4.000 C 
, 2^{ ‘._j.C 00 _C 


IFCmtuDA.GT .e.5)_Kkt:jjOAsh>FCOA-l.O 

JUM_UF ^C»>£^ CV FR SPLCT»4UM fclDTH 

KPs2.*>>l»TH*K FTSF 

'M'sKFnC-NP 
KMsKKTSFtNl 
IK*<l.LT.l) Mil 

K2 = ^FFC***P _ 

IKHi.CT.KPISF)' hJihPTSF 
pr»‘ Fso. 

“DC "2 CO" |sWl/0 
pr»th=p :*• H^i*SLT(i ) _ 

■(jC 2CI l = SM^ hPTSF - • - 

P=P':»»Kk«t<SLT( 1 ) _ 

PC»rs=P Jhtk/XNOkyOsf) 

■fXT~T HE' TH9TE~V‘ AVTAb‘LE'S'“l MC” A n”A F R A y 7~A NC‘ 
pf.AFS(i7J/isT)Vp‘n,rLP 

_ pf. r.(./Jy lCT)=FktUDA_ 

PE A> «;( j/j/IST)s*IDTH 
OUTPUT J/FPttnA 

HTUF *. - - -- 

f M 


SUfPCUTlNt PLSPfC ( 1ST) 


ORIGINAL -PAGc IS 
OF POOR QUALITY 


kOUTlNF TO PLOT Tilt SPtCTPUP 


PSLT 


hSLl AkkAY CONTAINING SPECTktk 

NPTS NUMhLR np POINTS IN PJLl 

1ST NUMbtk' OK ISTIPATOR 

IPLC NU'^bER OK KE.COPD 

LOGS'CL FLAG FUk LINEAR OR LOG 1LCT~ 
XNnPk ARRAY CONTAINING NOR AL I Z A T I CNS 


if i. 

1 j . 


I I 


Cll*'^CN/SIGNl/£IG(l024),kSL7(lC;i-n,W0PK(lC:4), 

iNPTSK,NPTS/NS1RT/NTUT/NJAT/PSIG/SNkF/SNh/ 

21 SYNC AT/GSIG/CNOISF/FREU/ SIG^ A _ _ _ 

C C ^ r N / P L S / X'N* C PV ( 1 6 ) / I C T ; L 0 G S C L , F I P'S 1 / fp X F L /L N A t ( 4 ) ^ L P L 

_CiP''CN NF ST/NAV/ IKSTy I ».NU*, LUk Dy ALPHAy kCTS( l2b)y NA 

COE ''ON ‘/£k/isk( o; 15)' ■ 

COkPLEX <;iGyWGTS _ 

LCCICAL FIRST/ IPL/PLCH/ IS* 
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K07.000. 
000 
H V.OOO 
It^O. 000 
1* 11. COO 

id ;.ooo 
id.3«ooo 

l^M. 000 






C 

C 


J)J‘'U.£10N . IPTO)^ IIH(3) . 
ISTTslST 

JK (lSh(l?JJ_JiQ10_95 

CbTU V7 

U (LLC5CL)~lCC^100/^00 


;f.t.C 00 .c_.„ 

ld.^.ooo 100 


LINEM) scale; 


)iJ7 * Cj2 

coo 


JT.b. 

Il V, 

i5:*c. 

;v:i. 

1 'L' ^ . 

1 L7 4 , y u 0 
JKt-.C00_ 
l#2t.COO C 


Y I N S 0 . 


ORIGINAL PAGE IS 
O F PQ Oft lAi iTv 


. 000 . 
i5:*c. coo 
coo 
coo 
x^n^Qo 
!l:4. 000 


V T I C s V A X / 1 C . 
5C TC 400 


LOG. SCALt 


200 


y»>^ysXMJKM( ISTT) 
yi- lN = -d0." 

_ Y'^AX^J. 


1^27. COO 

rT-f-«coo c 
1 J>y, 


Y *- A X = 3 . 

jnic = ic. 


00Q_C _ 
ib3C . coo C 
IlM.COO 300 
lJi2.C0) C 

1£33.C00 

1S34. COO 

1 ]3t.C00 

■]t3r.cdo C 


SfT Ur WINDOW 


_C«LL T.I'.DJ (‘'CylO00/?by 70C) 

NPTCE =/*M-7SF 

_ XMX = HtJAT(NnSF/2) 

Xk'i'jsi.-xyAX 

CA LL U . U DO ( X*^l Ny XVA X, YH I N^ If' A XJ 


JCO 


y\'ii,w.c 
1 |3e . coo c 
.1 J3^.c:o 
it<c. coo 
jl^i^coo. 

1 . coo 

15^3, coo 

1T4^. 000 
liV-. COO. 
Ib4t . COO 

iHt. 000 

it-;*-, coo 

IL50.C00 
ifit.COO 
coo 

iC* CuO_ 
000 C 

.lit V. coo c 
lb? f . coo c 


IF SWITCH 13 JS UP/ .00. AN AXIS. AM' NtSET IFXFL 


CQNTI'^U!- 

If ((.NUT.Fl'hST)) 
i.xi=i, 


CO TO 400 


» \] 


AX2iV-lN 

. AXJrFLUAl INPJSF)/IC, 

AX4 = YT 1C 
.. CALL AXlS(/.Xl/AX2yAX3^AX4) 


... ^ -I 


IPXFLsC 


-C- 

C 

c 


DtTfcN^'INL OUTPUT OF DASHED Llf.t 


400 


LslSTT-1 
I F_ X J b 7 T . F .J . P ) 
IF (Ibll.FU.O 
J f‘_U G U *i 2 ‘U. 


L = 3 
L = 4 
L = 0 


PLOT SP'CTpUM 


![ 


b7, CPU.. 
. COu 


M TS2 = bPTSF/2 _ 
y = hSLT( f.Pl 1 ) 


J 

' I 
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If-J^.coo 

litt.OOO 

1 fit . 000 
jl 5 ?.C 00 _ 
Jbtt.COO bOO 

IfrSS^COO 

\ 'l| 70 . c 6 o 

fl^^KCOO 

7 J^ 7 <,CwO 

i 1173.000 

1?7<.000 bOl 

147^000 

iPt.O'OO C 

1^77,000 


147^000 

iPt.O'OO C 

1^77,000 

lt)7b. COO 

JT7S.C00 C 

iJ^C.COO C 
l‘.«l.COO C 


. cco 

iJ-j.OOO c 

U 54 . coo C 

li-*‘. 000 _c 

j'J?t.COO c 

1*‘'7.COO C 

U5t.-000 

iJ»j 05.000 

live. coo 

\i M.COO 

iT' 7 ^coo 
: 4 ^ 3 . coo 

*:lye .coo 

i‘97.0C3 

li9t.C00 
1|QV.000 30 
■iTOl.COO 

ltd. coo __ 
'if oi.’coo c 

1*: 03. COO_C__ 

nOf.cc'O c 

J 'Ot.OOO 

.OtVCDO c 

1C07.000 C 

Y'otTcco c“ 

i os.coo_c 

"itioVcoo c 


X S 1 

Y»KS 

IKL 

CXLL 

CesT 

7s I 

Vs; 5 
iML 
CALL 
"CCM 
_OLL 
7 p 

KHU 

’lnD 


LT(NPTSFt_l) 

CCSCL.Cl.C) yslO.*ALO*ClO( V/XKkM) 

PASHA ^X/VyL) 

I SUl 

01 I s 1 ^KPTS 2 

I T( 1 ) „ 

Lrh'.L.'ll .0 VsiO.* ALUG 10 (Y />•.»« M) 

l AS-(A_^Xy Y_/IJ 

IMI 5 . 

A ». *■' 0 r> f 

T Lb s SH T SF72 

uu 


Sl'L'f-CUTl f.E 5PCAV (ISTT/ILPIO) 

"b L gcV~a' V f. P A c I ~K Tk s tHuTTs T rp'A* tTs 

b »f usb NTI AL_AVE»'ACE AFT^ ^F£T_^Ay. 
*HT IS HSLT FCR HLOITING 


CLHH 

JKFTS 

^ISYS 

CL 

F*.tV 

iVcT 
_XD=1 
XA='l 
_DC 3 
taCKK 
RSLT 
Fb*TU 
LSD 


cs/s 

DAT/ 

cs N 
MHLL 
= b LO 

rp 10 

. /FS 
.-XD 
0 1 = 
(l) = 
(li- 
ps 


ICNL/ S1C(1024)/PSL1( 1024)/SOKK(1C24)/ 

■»‘Sy ASJIPT/STbTyNDAT/PSIC/S,NKF^5NPy _ 

r,SIC/GsniSL/FHFi/SlC“C 

::STy» AVylESTy IhSUSyLCPr)yALP_HAy wCTS(l26)yKQ 
X SIC/WGTS 

AT_(SAV) 

.Lf.SAV) FNAV=FLuAT( ILFIO) 

AV 


lyNPTSF 
XA*WOPK (I) 
s-QPY ( I ) 


XD*PSLT(1) 


SUrPCI'TlSF CSKRIMI ST/ IHl ST) 
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I 


J.COO c _ 
13.000 C 
.•M3.000 C 

i M.COO C 
Ib.COO 
IC.OQO C 

i Jl.COO C 
IH.COO 
IS. 000 
.l 3 C • CuO 

r 21. 000 

23. coo 
.t?J. OCO 

1 34, COO C 
3t.COO C 
it2t,C00 C 

1 37^000 

ve.coo 
3S.CC0 100 

J n.ooo 
n.ooo 
J3.CC0 C 
U3j.CC0_C 
1 13 4 . C D 0 C “ 
li^t.CCO 
lt3t.C00 C 
i|37,c:d c 
ll3t'.c.<0 C 
LOS^cco.c 

1 40.000 C 
41. coo C 
43. coo C 

1443 . cco 

l|44.C00 

104' .coo 

iMf. OCO 

1 47,000 
4e. COO 

U4V.CU0 

It-'O. coo 

li *>1.1.000 

ltt'3 • OuO 

t SO.COO _ 

t4.C00 
‘-^r.COO 

t *)t . 000 
S7.C00 

, , ooo To 

:c‘>v. 000 
irt'».c4o c 

J.09. 100 
tt* 9.400 




KCOTINL tn ht/ib PEkKS,i'lklfS,ASD VAl-lANCES FFCH Fk;AKS 


CC“MCN f.t ST,1 AV / IFST/ I WSUnACRD/ALPHA/kCTSC I3fc)/N0 
Cd^-OS/HKi/PE AKSO/ 1024/ ») 

CU^'»'U^/:>IC^L/ S1C( 103^)/hSLT(1 034 )/ ( 1034 )/ 

_J ^PTS^ / KPTS^»SlPT/NTQT/f«UATyPSIC/SNHF/SM*/ 

2iSYM Al/GSTC/CSUISETFHFj'/SlCPA 

CDi'i'LJ * SlC^kCTS _ . _ . 

STOPL RUl'F.STLU .ISKO IN^OWK 

_LC_1C0 I = l/F»£T 

hrPM 1 )sirA'KS(IHrST/I/fST) 

CLMlN'Jt 

TURN 
1 NU 


_ S'PKUin.*: <TAT9 ( £ A N/ S T LL V / N h 1 N S ) 

. pnUlINi TO on STATIsnCAt analyses LN •CPK (MST) 

.nuTfUl j:_f Rn_'' 5T4 Ts 

XPtAN--> MEAN OF » 

SU'tV--> standard PEVIATION of k 


WORN 

F/SN 


(1C34 

R/ 


U ^ • V 4 

t ‘-3. Cl 
t4.C( 
‘-^r.Cl 


CCMPLLX_£ l.Cy hOTJ. 

SUMsO 

.. 5UPSD = 0 

DNst . 2H3l<j/F LOAT(NPINS) 

5U*'Cs0j 

SU^SsO. 

rC 10 I = 1/NF1N.S 

WrRSLK 1 ) 

X^rLi^-FLOAl (1-1 ) 

SU^C=SUMC4R*C0S( XN) 

‘•uy5s<;uk'n4>.«siN(Xh; 

CCNTINUF 

AN = ATAN2(SUHS/SUWC) _ 

OUTPUT XPFAN 
IF (XvkAN.CT.O) CniQ 91 
X*'tAS = -l*X*‘kAN 


hA,hCTS(l3e)/N0 



H 

u^o.oon 

lj-t)OV400 'c 

i' *>i. coo 9_l 
'lU2. 000*" 

cco 

"j| t iTtOO 

ltt4.0C0 

’ 16 ft. coo 97 
lj^t .ooo c 

■ U7C. 000 *98 
1671 , 000 

"ir7/*. COO c 

l|7j.CD0 _ 
U7< . COO* ■ 
It7‘.C00 _ 

ij 7f , COO* 
J177.000 
U7e .*C00 
lj7V.OOO 19 
IT eC. COO *■ 


i 


!’'• 

4 ► 

> M . 

I'** 

Jet . 

1 67 , 

I'"* 

t90. 
J^9l . 

:| 9 ;. 

o 94 . 
:T.9t. 
A 96 ; 
6V7. 
.? 9fc. 
i 99 . 
700. 

MO/. 
“7 03. 

} 704. 

;c‘ . 
vot . 


COO 
000 
COO c 


cco 

coo 

coo 

CGO 

000 

000 

ooo' 




1707. 


f; oVi 

i_j40S. 
171C. 

III. n'.'- 

l.H /. 


c 

c __ 

c _ 
c* 

c 

*c 
coo c 
coo* c 

coo 

000 

coo 

000 

000 

coo 

coo 

coo To' 

coo 

0*00 

coo c_ 
000 c 
coo 200 
coo' ■' 

ooo_c 

cuo'c 
000 c 
coo oo'o* 

000 

coo ■ 


_X»'fA\*6.:i>3l9-XH£AN_ 
OUTPUT X^'KAK"** 

DU 97 1*1/NHNS 

XW*l.««KLCAT( r-l ) 

Jl^r »K*'t 

ir(x*-D.rT.o*)iNL»Exri-i 
IMAG'D. LT.O)COTO__^ 
CONTIKUI 

__ CUTPUT_ IXDr X^OW 

|HDEX*lNCEX-9blKS/2-l 


^’^Al p7. ■ ■■ 


XvrAf.s(lMiLX-.5)/FLUAT(SH5NS) 
CUTPUT'XPf ND'EX 

LO 19 I=l,Nt$T 

*T£STs»'J^K (1 )-XHF.AS 


sATf Sl**2 

US7 = Su9T( AUST) 

USTlsl .-nsT 

1 r ( T t S T ."G T*.*.T ) A“T E ST^f t*S T I • • 2 
SU'^St sATESTtEUMSD__ 

"* STCFWsSwPl (SUMS0/(NEST-1 ) ) 

hfcTURN _ 

*£f.D 


SUhhCOllHE *kl*lST"(HVAX;96lNS/ IHIST)' 

KUU 1 1M~TU~"UU* HI STOCK A9. ON WCKK(f.FST) 

|iCf-K--^NPUT ARRAY *1TH PAN [ATA 

K£LT--> UUTP’UT* ARRAY ■|TiTH*"H1S1CCRA> 

lajcest value. I n. hi ST 

IGNL/ S1C( 1 024),PSLT( 1024),WUPK( 1C24)/ • 
lMTSF/hPTS/NST«T/NTUT/NOAT/PSIC/SNPf /SNR/ 

21 rVNrAT/CSIG^£NClSE./mG/SJ[G^(. „ _ _ 

■ *C C R V C N N E'S T7 N A*V*/Ye‘S T / l <» N D W / L O R D / a L P*H"a7 r c T S ( I* 2*fc*) / n'o* 

CU»'I-LEX_S1G^NCTS 

■RMyr]. 

DC 10 I = l/NfclNS 

RSLT(I)sO. 

_DL 100 IRONIST 

IF "( lrilST-2) 200/ 300/400 

*r(NFR NJRVA*lTZE*D~TO**To*72.0) 

INDh y = 0.5* KU: AT (NM NS)* *nRK( 0.5 

GC TC 500 ■ 

"f f r‘7uV*. c y“nT r*^alized fCT (o~~ 5T*~c“'5 ) 

■*M.Dr‘XsFLCJAf(NtIN"S>*NORK*( I )^oT5* 

r.DF) = INULX ♦_Nb_lNS/^_* 1 

* CO Tr 500 


• — -1 


-I 



J 




I 

I 




(I 

i 

, ! 

1 1 
>1 


J 


i 

■ 

^ 1 




1 


I 


% 
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1 


.nu.ooo c 

! *»m,coo c 

'lb. coy c 

/It. 000 400 
17 17. c o o bOO 
T'lt-.OOO 

I/I9.CO0 

172C.C00 

{ '21 .000 

' 22.000 100 
_l7 <:ii 00 O.£ 

f »24.C00 

»? b . 000 

l72t.C00 C 

WT 27, COO.C 

1/2*-. CCO C 

_ r/:*s.oco 

’L7 3 bVoo 6 c 

1/31.000 C •*_••• 
//31.C00 C 

_1733,C00.C 

:134.C00 C 

l73t .ecu C 


fclPTH NnP*»*L121.0 TO (O.O/O.b) 






1SDEX*:.*FLCM (NhIKS)*«JM((I)40.5 

U ( I M • K . Li. 1 ) lS\>f tT\ 

IMl /*f'x X.CT.Nl'lNS) lNLt)C*NSINS 

fSLT( lMit*)«KSLT(lVDtX)„* 1 _ 

IK/'SLTC 1M)M).CT.H*«AX) HV * XsR SLT ( I NPF X ) 
IF (KSLT(UDEX2 tF(,,H.U.X) JlNDsISDEX. _ 
CevT INUI 

_ C U 7 t U1_H ► X x^ HIV’U 

1 - / T'j F N 


_f Ui^CUT lr.E hFLL T feK/STrSV^NblNSy I H 1 ST^l STJ,. 


kUJllNE TO PLOT hISTOCPAMS FPCm kSLT 


:f37.cco c 
:J»3t.OjO C 


, j 73 t . COO_C 

I ’'’4C.C00 

' iKpOO . 

1/42. CCO 

. v^v.ooo 

|/44 . COO 

. A 41.000 

174< . CuO 

l'47iC00 

i ~ 4 P.COO 

174S.COO 

.f bO.CCO 

..j'^:i.coo 

m 2 , oc-o 



:f -b4 . coo 

r / bt.COO 

L7bt . COO 
] i '*- 7 . CC 0 

I ^ VSt.OUO 
_J7‘ S. CCO 

C ' ' C . 000 

' cl^COO 

l7t.’. 000 lu 

i *'* ‘.COO 

't 4 .000 


CCHt'CN NtST/XAV / JEST/ U*JL)»/LOPD/ ALPHA,*C1S( 12E)/N0 

^CCy^'CSySJCNLy _SlCil0i4Jy PSLTC»Ci4)y kpPK(Ui4J^ 

lAf TSF / f.Pl S/ A‘lPT/ML-T/SOAT/PSlCy SNHF/SNP/ 

2lSV^DATy CSIG/CNOI5t/FRL3/SlCPA . _ 

CCVKCS/PLS/ XNOPM( 10)/ ICT/LOCSCL/FIPST/IPXFL 

_CC^‘'CN/LPLS/LAP.tL(10/a) 

CC^PLFX SIC/WCTS 

l;l-KNSION lCT(i)^ — 

IPC=(H^'AX/rLCAT(KF.!‘>T))*r0 0.^l. 

PC = i- LL'AT( IPC)*J-LOAT(NKST)/lOO. 

CALL (1./FLOAT(N'11NS)/0./PC) 

__CALL T^INL'O. ( 100/ 1000/ 100/750) 

PCl(; = f'C/10. 

P.7 SI C5»- LiLATLNJJ ^5JyiiL• 

XAXsl . 

_ IF (1H1S1.EC.2) XAX=NBi/S/2U 

CALL AXIS (XAX/0. /PTSIO/PCIO) 

CALL .PF.CUVR . .. . 

CALL X'UVEA (l./b.) 

_r.L 10 i-i/'iiliNs 

CALL LKAnA (F LUAT( 1 )/PSLT( I )) 

_CiLL T/F A4 A. (HUAT (!♦!)/ •'SLT(l)) 

CLNT I f.UF 

..CALL aevt 

CALL Afi-iJUt 


L 


I7tt.000 X 


j Yt.oo. 

11 7,001 

1 ‘Ab’. oo( 

if 0 . 

ini, 

iT 3 . 

11 4 . 

1 _ 77 *^ 

1 Vt .( 
V 17 .J 
I77t. 

1 VS. 

ii.c. 

itVk 

l^i. 

n*c 4 . 

j'T t.cc 
J 7.c: 


) no 


ACCtrT_'<7>*^lPU,H 

fcMTk '(lOVlCSO) LOHD/NbUS 


000 

ooo“ 

jVSjjCOO 
11 0 . 0 ^ 0 “ 
_ coo 

177;.'C00 
000 
coo 

ilU.CCO 

cco" 
_000 
I77t. coo 

CC-0 
ll.C.COO 
000 
o'oo' 
coo 
000 
coo 
coo 
_ :-o 

s. cco 

C.CuO 

1 ._coo 

nooo 


1090_ KCb**!! (' (CbDLB«*^ 13 

C ■ "cc TO ioo 

C13 C »M T fc XI O^ic 3_CJ IpPp ^ tlPt U^Kr IK S 


'J 'J 




1030 FCK^AT (' (CbUbs%13/* ALP*'/F4.2/' P1NS*'/H^')') 

JOO WMTF (10/2CCC> _NtSTyKTOT/SS^XbEAN/STDEV 

"2000 FCb^'lT (' ESTIMATES */I4/' SAMPLl S/E ST I MATE '/F5.1/ 

l_* Oil HEAN'n^Fb.5, * ST. DE V. : F 7^ 5^ Z) 

duTPOT • ' 

OUTPUT • • 


1F( IST.tU.l )cu:put 
1F( IST.' W.t lOUTPUT 
1F( I'JT.L'./.'t )OUTPUy 
IF( ! *;T.h .,t )CU:^UT 
LUTP'IT ' * 
OUTPUT ' * 
LUTPUT'HAkC 
INPUT «OAN 


'FfT HISTCGKAM' 

'hGTS PISTOCKA-' 

■•alpf'histccka«' 

'•'IITILAC HSTOCRAV' 


IKHCAN.EC.DCALL 

CALL ‘,'*PAC 

P'F TUPS 
I SO 


copy l-VES/O-Nu* 

HDCJPY 


C 

C 

C 


SUPPCUTiNt PFT 


1 “. 


iTs^V 

iV'* 

rl<t . 
rj*y 7 _. 
17 V>-. 

1 l-A* 

1 4jC. 
I h 0 1 . 

1 lii* 

1 1 ) 3 . 

Ib04. 
1 V-S 

iIdc. 

1^07. 

'I 


ipoe. 
1 j7S. 
10 . 

ilLL*. 

2 . 


ooo_c__ 

OOO C 

ooo_c 
coo c 
coo_c 
coo c 
0 y 0 _c_ 
coo c 
coo_c_ 
coo c 
ooo.c 
coo c 
coo c 
000 

O'f'J 

coo 

Ono_ 

00 ) 

OGO_C 

cob ci 


3. coo 

4. coo 

‘._000 

t.coo c 


UlSCPlPTlbN CF PARAMETERS: 

UC-AS INHUJ _C0NTA1 KS_COMPLEX _ONE-^lMiNnCJAL 

AKKAV TU bF. TRANSFORMED 

-AS 0U7PUT-CCNTA1SS_THE COMPLEX FUUPJJ,P JPjkNSrCbb 
M-DH-INPS TPE DIMENSION CF A. 

IliL JJ ^ CjyF_AJ M 

N-L’l'yh'sS'lON CF A 

TNVr.rPllUN FAKAMFTEk 

1. 0 'TA“kE S TR ANSFURM~CF A 

.-1 I) JAXE§.J.NVi:RSi TRANSFOKM _0F A 

NOTE: N MUST bE A F’UKEK CF 2 


CCMhCN/SlGNL/ S 1C( 1 024 ), RSLT ( 1024)/ iORKC 1024)/ 

INPTSF/NPTSy KSTRTyNTiJTy N J A Ty PS I G/ Sf^W 

2*ISyMiAT/GSIC,CNUISE/FRE”j/ sigma 

__C.C P L E X SI Gy Uy Ky 

N = i! • • M 

lie 2 I H)=l /N 

nufpuT sicc'nc) 

NV2 = N/2 

NM = r .- 1 

_Jsl 

OUTPUT N/NV2/J/T.NV 



134 


THIS LOUP THt 

.CPtHATiON UK THf iNUtX Ur A . 


PC. 7 - 

IH ( 1 .Gc.* J) CO TO b 

_ TrSIC(J) 

S 1 C( J)sSIG( 1 ) 

_ SIC(I)sT . 




.. • ^ • - JJL 


103 .C 0 O 

Vb* 3 A.'t 00 ' 


byit‘ 3 t )9 

,w_ 

1 /M 


•* 93 t . 000 

1 H 37 , coo 
IbBe.COO 
ife 39 jC 03 . 
1 - 4 C. CCO 


CuScVl/Ltl )7 tkV* SI N(P 1 /Lt 1 ) ) 

■■ijl.U 

■ J / N / L L 


;( iF)*u 

P)=5IG(J)-T 

)-ilG( 1 )♦! 

» _ _ — . 

; - 1 N V 

JUTPUT K^KV2yJ/T^y 
»U /fl IC 1D=1/N 

jUTPUT i JGtlDlO) 

NV .GT. C) RHUftK 


IbSS . COO 
I cS 4 . CjO 


NVi -1 


HAVE THE INVEKSE TRANSKOf-M 


S=l . /U 

C=C^PLX(?, 0,3 

UC 30 1 =1/ N 

I ) = C.*S1C(.I ) 
OUTPUT SIC(I) 
COKTINUl 
pi TURN 

i r.lj 


__ f.UMHOUTlNt AX 1 S(X/ Y/XTICy YTIC) 
CALL SfFP^OPlNyX.yAXyYwiNyVPAX) 

CALL VbVEA(XMlNy Y) 

"call LI- Ai^A(XVAX,Y) 
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t 


•« • 


. C) CO TO^ 
IK)/XTIC^1 
>/>T 1C 
• X* 

( 

>C^V) 

[O/lO) 
[0/-20) 


10 





• .LE . 0) CO T.l \0 
•XX-NM1‘. )/VTlC^ 1 
''U)/VT1C __ 
C'HKY 

1/NV_ 

LA (X/VO) 

«EL (10/0) 

.wlL (-;0/0) 

1C 


nf'N ii<nn(j/NU) 


lHTI-slhlTl.*2MJl-^*J2> 
2C0 JlsJ^ 

1. f T I K f. 


EUKCTION LKANCDUM) 


K ANDOH 


« h£TfcFtN 0 ♦ I 


CM»*1 NSinX UAN(1 ) 

NS) 1 Dsl 

CALL kSUU(MK/ NSt'tP) “ 
r. S S 2 1 0 0 • K A N ( 1 ) 

TALI ivM'Ua AN/S*-! KD) 

L' C 2 - 1 • * ( N S S • >< A N ( 1 ) ) 

’call )*iUUO.AN/<*SLtU) 

D ^ A L S • »- A N ( 1 ) 

KETUKN 
\ NP 
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